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ABSTRACT 


Linearized  six  degree  of  freedom  rigid  body 
aircraft  equations  of  motion  are  presented  in  a 
stability  axes  system. 

Values  of  stability  derivatives  are  estimated 
for  two  representative  STOL  aircraft  -  the  DeHavilland 
of  Canada  "Buffalo”  and  "Twin  Otter".  These  estimates 
are  based  on  analytical  expressions  included  in  the 
report.  The  combination  of  the  equations  of  motion 
and  the  estimated  stability  derivatives  provides  an 
aircraft  model  which  is  useful  for  Navigation,  Guidance 
and  ATC  Studies . 

Resulting  transient  responses  to  control  inputs 
are  presented. 
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SYMBOLS 


i()  j()k() 


Ix,Iy ,zz 


Airplane  lift  curve  sloDe 


rad 


lift  curve  slope  of  surface  () 

(of  wing  when  no  subscript) 

2 

aspect  ratio  of  surface  ()  =  b,*/S() 

(of  wing  when  no  subscript) 

span  of  surface  ()  (of  wing  when 
no  subscript) 

mean  chord  of  surface  ()  (of  wing 
when  no  subscript) 

nondimensional  stability  derivative 
(defined  in  Section  VI) 

parasite  drag  coefficient  of  aircraft 

lift  curve  slope  of  wing 

aircraft  parasite  drag 

wing  efficiency  factor 

force  component  along  the  ()  axis 

gravitational  constant  =  32.2 

altitude 

height  of  fuselage  at  wing  root 

CG  position,  fraction  of  c 

neutral  point  of  aircraft,  fraction  of  c 

component  of  angular  momentum  along 
the  ()  axis 

unit  vectors  along  the  X,Y,  and  Z 
axes  of  the  ()  coordinate  frame, 
respectively 

aircraft  rolling,  pitching,  and  yawing 
moment  of  inertia,  respectively  (further 
defined  in  equations  13-15) 


rad 

ft 

ft 


rad 

lbs 


lbs 

2 

ft/sec 

ft 

ft 


slug-ft  /sec 


slug-f t 


2 


Jxy 

product  of  inertia  =  /xydm 

slug-ft 

Jyz 

product  of  inertia  =  /yzdm 

slug-ft 

J  =  J 

zx  xz 

product  of  inertia  =  /xzdm 

slug-ft 

L,M,N, 

scalar  components  of  the 

applied  external  moment  along 

the  XA,  Y  ,  and  ZA  axes,  respectively 

ft-lbs 

l() 

di stance,  quarter  chord  of  aircraft 
mac  to  quarter  chord  of  surface  () 

ft 

m 

mass  of  aircraft 

slugs 

P,Q,R 

scalar  components  of  the  angular 
rotation  vector  of  the  aircraft 
along  the  XA,  Y  ,  and  Z  axes, 
respectively  A 

rad/sec 

AP  ,  AQ, AR 

perturbed  portion  of  P,  Q,  R, 
e.g.  AP  =  P  -  PQ 

rad/sec 

q 

dynamic  pressure  =  ^  P 

lbs/ft2 

s 

J 

laplace  transform  variable  =  ^ 

sec“l 

s() 

area  of  surface  ()  (of  wing  when 
no  subscript) 

ft2 

t 

time 

sec 

AT 

change  in  thrust  due  to  pilot 
throttle  input 

lbs 

U,V,W 

scalar  components  of  velocity  of 
the  aircraft  along  the  X  ,  YA  and  Z 
axes,  respectively  A 

ft/sec 

AU,AV, AW 

perturbed  portion  of  U,  V,  W, 
e.g.  AU  =  U  -  UQ 

ft/sec 

Uo 

equilibrium  or  reference  value  of  U 

ft/sec 

VR 

resultant .velocity  vector  of  aircraft 
=\/u2  +  v2  +  w2 

ft/sec 

W 

width  of  fuselage  at  wing  root 

ft 

w 

weight  of  aircraft  =  mg 

lbs 

X,Y,Z 

scalar  components  of  the  applied 
non-gravitational  external  forces 
along  the  XA,  Yft,  and  ZA  axes, 
respectively 

lbs 

IV 


v  y  Z  axes  defining  the  ()  coordinate 

0  0  ^  frame 


AX  AY  AZ  perturbation  portion  of  X,  Y,  Z, 

'  '  e.g.,  AX  =  X  -  XQ 

x,  v,  z,  distances  along  the  XA,  YA,  ZA 

axes,,  clarified  by  subscript 

y-^  distance,  aircraft  centerline  to 

inboard  end  of  aileron 


a 


Y 


distance,  aircraft  centerline  to 
outboard  end  of  aileron 

angle  of  attack  =  tan  ^  W/U 

angle  of  sideslip  =  tan  ^  V/U 

deflection  of  control  surface  () 
(positive  deflection  produces  positive 
moment) 

perturbation  portion  of  6,., 
e.g.,  A6()  =  6()  -  6()^ 

flight  path  angle  of  aircraft 
=  0  -  a 


r 


£ 


n() 


0,$,4' 


A0,A*,AY 


P 

a 

T 

Subscripts 

A 

C 


dihedral  angle 
downwash  angle 

efficiency  of  tail  surface  () 

Euler  angles,  defined  in  Figure  2 

perturbed  Dortion  of  0,  $ ,  f , 
e . g . ,  A0  =  G  -  0^ 

atmospheric  air  density 

sidewash  angle 

aileron  effectiveness 

aircraft  body  coordinate  frame 
Earth-aircraft  control  coordinate  frame 


E 


Earth-centered  coordinate  frame 


I 


inertial  coordinate  frame 


L 


Earth  local-vertical  coordinate  frame 


lbs 

ft 

ft 

ft 

rad 

rad 

rad 

rad 

rad 

deg,  rad 
rad 

rad 

rad 

3 

slug/ft 

rad 


v 


flaps 

fin 

due  to  gravity 
stick-fixed  neutral  point 

horizontal  tail 

equilibrium  or  reference  condition 


wing 


1.0  Introduction 


There  is  no  such  thing  as  an  exact  mathematical  model  of 
any  physical  phenomenon.  All  mathematical  models  are,  therefore, 
approximations  to  reality;  approximations  based  on  an  assumption 
set  which  should  be  clearly  recognized  by  both  the  writer  and 
the  user  of  the  equations  comprising  the  model. 

The  basic  set  of  assumptions  must  always  be  derived  from  the 
desired  application  of  the  model.  To  develop  and  utilize  a  more 
exact  model  than  that  required  for  the  job  at  hand  is  to  pay  an 
unnecessarily  high  price  in  man-hours  and  computer  time  both  in 
the  use  of  the  model  and  in  the  gathering  and  formatting  of  the 
data  required  by  the  model. 

The  model  presented  herein  is  a  linear  perturbation  model. 

It  was  developed  for  use  in  exercising  sets  of  4-D  guidance 
equations  which  are  being  developed  for  application  in  STOL 
terminal  area  guidance.  The  aircraft  which  are  modeled  (the 
C-8  and  the  Twin  Otter)  were  selected  as  representative  of 
two  classes  of  aircraft,  i.e.  light  and  medium  STOL  propeller 
aircraft,  and  are  of  interest  only  as  representations  of  the 
classes  from  a  guidance  and  aTC  viewpoint. 

This  report  is  submitted  as  partial  documentation  of  work 
carried  out  in  support  of  PPA  18-0,  dated  December  1,  1970. 

In  particular,  it  documents  the  work  performed  under  Task  2 
of  this  PPA  for  the  "Buffalo"  and  "Twin  Otter"  aircraft. 

(Similar  data  are  required  under  this  task  for  a  third  STOL 
at  a  later  date.) 

The  applicability  of  the  present  model  is  discussed  in 
detail  in  the  next  section,  Section  2.0. 
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The  next  3  sections  are  devoted  to  the  derivation  of  the 
required  linearized  equations:  In  Section  3.0,  required 
coordinate  frames  are  defined.  Section  4.0  derives  the 
generalized  kinematic  equations  of  motion,  utilizing  very 
few  approximations.  Finally,  in  Section  5.0,  the  desired 
linearized  equations  are  developed  from  the  generalized  ones. 

Section  6.0  presents  analytical  expressions  for  the 
required  stability  derivatives.  These  expressions  are  used 
to  generate  numerical  values,  given  in  Section  7.0,  for 
two  representative  STOL  aircraft. 

Supporting  material  is  included  in  the  Appendices.  In 
Appendix  A,  analytical  expressions  for  selected  stability 
derivatives  are  developed.  Appendix  B  contains  estimates  of 
moments  and  products  of  inertia  for  the  two  aircraft  under 
consideration.  Appendix  C  contains  a  step-by-step  calculation 
of  the  stability  derivatives  summarized  in  Section  VII. 
Appendix  D  presents  typical  transient  responses  to  control 
inputs,  calculated  using  the  equations  and  derivatives 
developed  in  this  report. 
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2.0  Applicability  of  Mathematical  Model 

As  stated  in  the  Introduction,  a  linearized  representation 
of  aircraft  motions  is  required. 

The  model  is  intended  for  use  as  a  tool  in  the  preliminary 
design  and  analysis  of  STOL  aircraft  control,  guidance,  and 
navigation  systems.  In  this  phase  of  analysis  and  design,  the 
unstabilized  response  of  the  vehicle  is  adequately  established 
by  means  of  a  linearized  analysis.  The  linearized  model  lends 
itself  to  the  use  of  such  techniques  as  root  locus  analysis 
and  frequency  domain  analysis  . 

Many  approximations  are  required  to  develop  a  linearized 
set  of  equations  from  the  generalized  kinematic  equations  of 
motion.  All  assumptions  used  in  the  derivations  are  explicitly 
stated  in  Sections  4.0  and  5.0.  They  are  introduced  as  they  are 
needed  and  are  consecutively  numbered  so  that  the  reader  can 
easily  establish  the  degree  of  simplification  at  any  point. 

The  major  simplification  is  the  introduction  of  small 
disturbance  (or  "perturbation")  assumptions.  Under  these 
assumptions,  aircraft  motions  are  restricted  to  small 
excursions  -  perturbations  -  from  an  equilibrium  flight  condition. 
The  major  virtue  of  this  assumption  is  that  it  vastly  simplifies 
the  equations. 

Its  use,  of  course,  limits  the  applicability  of  the  equations 
to  a  certain  extent.  The  reader  is  cautioned,  therefore,  to 
determine  the  effect  of  this  assumption  (and  of  the  others)  before 
applying  the  equations. 
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A  reservation  must  also  be  stated  concerning  the  stability 
derivatives  presented  in  Sections  6.0  and  7.0.  These  derivatives 

i 

are  not  based  on  wind  tunnel  or  flight  tests,  because  data  from 

,  ,  t  i  I  . 

these  sources  were  not  available.  They  have  been  developed  by 

'  '  i  '  . 

analytical  niethods ,  and  have  been  augmented  in  some  cases  by 
generally-applicahle  empirical  data. 

,  I 

•  '  * 

The  intended  purpose  of  the  estimated  derivatives  is 

i 

to  establish  in  representative  fashion  the  dynamic^  of  small 
and  medium  types  of  STOL  aircraft.  They  should  not  be  used 

'  ’*  •  l 

••  'I 

as  the  basis  for  an  evaluation  of  the  flying  qualitites  of 

•  i  * 

the  "Buffalo”  or  "Twin  Otter"  or  qf  the  suitability  of  these 
aircraft  for  any  specific  mission.  1 

I  1 

•'  i  ;  '  .  '  1 

1  i 

■  i 

I 
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3.0  Definition  of  Reference  Coordinate  Frames 

The  equations  of  motion  of  an  aircraft  are  based  on  Newton's 

I  l1  f  i  I 

second  law.  This  law  relates  the  forces  applied  to  the  aircraft 
to  the  acceleration  (or  change  of  momentum)  of  the  aircraft  with 
respect  to  inertial  space.  It  is  usually  convenient  to  define 
applied  forces  and  moments  with  respect  to  a  frame  fixed  in  the 

aircraft.  Further,  aerodynamic  forces  depend  on  the  motion  of 

1 

the  aircraft  with  respect  to  the  air  mass.  Finally,  the  motion 

of  the  aircraft  with  respect  to  the  Earth  is  frequently  of 

/ 

interest.  Thus  it  can  be  readily  seen  that  several  reference 
coordinate  frames  are  required  to  completely  describe  the 
applied  forces  and  the  resulting  motions  of  the  aircraft. 

Reference  coordinate  frames  to  be  used  in  this  analysis  are 

J  i  I  i  1 

defined  in  this  section.  Insofar  as  possible,  axis  systems  have 
been  defined  so  that  senses  of  rotation  and  translation  are 

i 

similar  for  small  rotations.  Positive  force,  moment,  and 
motion  vector  components  are  ddfined  to  be  in  the  positive  sense 
of  the  axis.  To  the  largest  extent  possible,  the  symbols  and 
conventions  used  are  consistent  with  those  in  common  usage  in  the 
guidance  and  control  fields  and  with  those  used  by  NASA  for 
aircraft  stability  and  control  work. 

The  Inertial  Coordinate  Frame  (I)  will  be  defined  first. 

This  frame  is  nonrotating  with  respect  to  inertial  space.  The 
origin  is  the  center  of  the  Earth,  with  the  Z-j-  axis  coincident 
with  the  Earth's  axis  of  rotation.  The  and  Yj  axes  then 
lie  in  the  equatorial  plane.  It  is  assumed  that  the  linear 

'  i 

and  angular  accelerations  of  the  Earth  with  respect  to  inertial 
space  as  it  moves  in  its  solar  orbit  are  not  of  interest.  This 
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coordinate  frame  (as  well  as  the  E  and  L  frames  following) 
is  shown  in  the  sketch  of  Figure  1. 

The  Earth-Centered.  Coordinate  Frame  (E)  is  fixed  with 
respect  to  the  Earth.  Its  origin  is  at  the  Earth's  center  with 
the  ZE  axis  coincident  with  the  Zj  axis.  The  X£  and  Y£  axes 
lie  in  the  equatorial  plane,  intersecting  the  Earth's  surface 
at  convenient  points.  The  E-frame  can  be  chosen  to  coincide 
with  the  I-frame  at  a  particular  instant  of  time. 

The  Earth  Local-Vertical  Frame  (L)  is  a  local  geographic 
frame.  Its  origin  is  the  center  of  mass  of  the  aircraft  with  Z 

Li 

along  the  vertical  defined  by  the  local  gravity  vector  (positive 

downward) ,  X  parallel  to  geographic  North  (positive  to  the  North) , 
L 

and  parallel  to  geographic  East  (positive  to  the  East) . 

The  Aircraft  Body  Coordinate  Frame  (A)  is  fixed  to  the 
aircraft  and  rotates  and  translates  with  the  aircraft.  Its 
origin  is  the  center  of  mass  of  the  aircraft.  The  XA  axis 
is  chosen  in  a  convenient  forward  direction  in  the  plane  of 
symmetry.  (The  exact  X  axis  location  is  specified  in  Section  V.) 

A 

The  Ya  axis  is  normal  to  the  aircraft's  plane  of  symmetry 
(positive  to  the  right) ,  and  the  ZA  axis  is  in  the  plane  of 
symmetry  (positive  downward)  and  orthogonal  to  the  XA  and  Ya 
axes.  The  A-frame  is  related  to  the  L-frame  (and  to  the  next- 
defined  C-frame)  in  Figure  2. 

The  Earth-Aircraft  Control  Coordinate  Frame  (C)  is  also 
centered  at  the  center  of  mass  of  the  aircraft.  The  Z q  axis 
is  aligned  with  the  local  gravity  vector  (positive  downward)  and 
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is  therefore  coincident  with  the  Z^  axis.  The  axis  is  the 
intersection  of  the  horizontal  XL-YL  plane  with  the  vertical 
plane  containing  the  axis.  The  Yc  axis  completes  the 
orthogonal  right-hand  system.  The  C-frame  is  an  intermediate 
frame  needed  to  define  the  Euler  angles  describing  the  relation¬ 
ship  between  the  Earth  local-vertical  <L)  frame  and  the  Aircraft 
body  (A)  frame.  In  their  order  of  rotation  (which  must  be 
preserved)  the  Euler  angles  are  defined  as: 

1.  Heading  (Y) s  angle  of  rotation  about  ZL 
from  XL  to  XqJ 

2.  Pitch  (0) :  angle  of  rotation  about  Yc 
from  Xc  to  XA; 

3.  Roll  (<t>)  :  angle  of  rotation  about  XA 
from  Yc  to  YA< 

These  Euler  angle  rotations  are  shown  in  Figure  2. 
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4.0  Derivation  of  Generalized  Kinematic  Equations  of  Motion 

The  generalized  equations  of  motion  are  obtained  (as,  for 
example,  in  References  1  and  2)  by  equating  forces  and  moments 
acting  on  the  aircraft  to  the  rates  of  change  of  linear  and 
angular  momentum  of  the  aircraft  with  respect  to  inertial  space: 

(1) 
I 


F  =  (mv) 
eft 


M  = 


jd 

dt 


(2) 


In  these  equations,  F  is  the  force  vector  and  M  is  the  moment 
vector  acting  on  the  aircraft.  Linear  and  angular  momentum  vectors 
are  represented  by  mV  and  H  respectively.  The  subscript  I 
indicates  that  the  time  rate  of  change  of  the  momentum  vectors 
is  with  respect  to  inertial  space. 

It  is  of  more  interest,  however,  to  express  these  momentum 
changes  in  terms  of  an  axis  system  that  is  fixed  in  the  aircraft 
and  that  is  therefore  translating  and  rotating  with  respect  to 
inertial  space.  This  axis  system  is  the  A-frame  defined  in 
Section  3.0.  Equations  1  and  2  can  be  expressed  in  A-frame 
coordinates  as 


F 


x  mV 


(3) 


M 


+  a)  x  H 


(4) 
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In  these  equations,  the  subscript  A  indicates  momentum  changes 
with  respect  to  the  rotating  A-frame.  The  cross  product  terms 
account  for  the  fact  that  the  A-frame  is  rotating  at  a  rate  w 
with  respect  to  inertial  space. 

Equations  3  and  4  are  completely  rigorous.  At  this  point, 
however,  it  is  convenient  to  make  several  assumptions  to  facilitate 
further  development  of  the  equations  : 

Assumption  1:  The  mass  of  the  aircraft  does  not 

change  significantly  in  the  interval 
of  interest,  that  is 


af  (mv)  ■  ”  at  (^> 

Assumption  2:  The  rotating  earth  can  be  considered 
an  inertial  frame  for  the  purposes  of 
this  analysis.  Therefore  the  E-frame 
is  assumed  to  be  an  inertial  frame. 

Assumpiton  3:  The  aircraft  is  a  rigid  body.  The 

contribution  to  R  of  spinning  propellers 
can  be  neglected.  Control  surface 
dynamics  need  not  be  considered. 

Assumption  4:  The  YA  axis  is  a  principal  axis  so 

that  the  products  of  inertia  JXy 

and  J  are  zero, 
y  ^ 


By  virtue  of  assumption  2 ,  the  vectors  V  and  5  of  equations 
3  and  4  are  the  motions  of  the  aircraft  with  respect  to  the  Earth. 
These  vectors  are  further  defined  in  terms  of  their  A-frame 
components . 


V  =  u 

aa 

+  V 

3A 

+  W 

(5) 

o. 

ii 

1  3 

rA 

+  Q 

^A 

+  R 

(6) 
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where  iA,  j  ,  and  are  unit  vectors  along  the  Xft,  Y A,  and  ZA 

axes,  respectively.  Similarly,  F,  M,  and  H  of  equations  3  and  4 
can  be  expressed  in  their  A-frame  components: 

F  =  Fx  *A  +  Fy  Ja  +  FZ  *A  (7) 

M  =  L  iA  +  M  jA  +  N  kA  (8) 

5  -  Hx  h  +  Hy  Ja  +  Hz  EA  <9i 

Development  of  expressions  for  the  scalar  components  of  H  is 
rather  lengthy  and  will  be  omitted  here.  Following  the 
derivations  in  Chapter  1  of  Reference  1  or  Chapter  4  of 
Reference  2,  for  example,  produces  these  expressions  for 


the  components  of  H : 

Hx  =  PIx  -  QJxy  -  RJxZ  (10) 

Hy  =  QIy  "  "  PJxy  (ID 

Hz  =  RIZ  -  PJXZ  -  QJyz  (12) 

where  the  moments  and  products  of  inertia  are: 

Ix  =  '(y2  +  z2)  ^  (13) 

Iy  =  / ( z2  +  x2)  dm  (14) 

I  =  /(x2  +  y2)  dm  (15) 

z 

JXy=  /  xy  dm  (16) 

JyZ=  /  yi  dm  (17) 

Jxz=  /  xz  dm  (18) 


and  x,  y,  z,  and  the  mass  element  dm  are  defined  in  the  sketch: 
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? 

z 


By  virtue  of  Assumption  4,  equations  10,  11,  and  12  reduce  to 


PIX 

■  RJXZ 

(19) 

QIy 

(20) 

RIz 

-  PJXZ 

(21) 

Now,  by  using  equations  5  through  9  and  equations  19  through  21, 
equations  3  and  4  can  be  expanded  to  give  A-frame  components  of 


aircraft  accelerations  with  respect  to  the  Earth: 

Force  Equations 

Fx  =  m  IU  +  QW  -  RV]  (22) 

Fy  =  m  [V  +  RU  -  PW]  (23) 

F  =  m  [W  +  PV  -  QU]  (24) 

z 

Moment  Equations 

L  -  Ix  P  +  (Iz  -  Iy)  QR  -  Jxz  (R  +  PQ)  (25) 

M  =  I y  $  +  (Ix  -  Iz)  RP  -  Jxz  (R2  -  P2)  (26) 

N  =  Iz  R  +  (Iy  -  Ix)  PQ  -  Jxz  (P  -  QR)  (27) 
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The  A-frame  forces  (Fx,Fy,Fz)  and  moments  (L»M,N)  in  equations 
22  through  27  represent  all  of  the  external  forces  and  moments 
acting  upon  the  aircraft.  These  forces  and  moments  are  due 
to  aerodynamic  loads,  control  and  propulsion  systems,  and 
gravity . 

Gravitation  Forces  -  The  gravity  force  is  a  vector  quantity 
of  magnitude  mg  acting  along  the  positive  ZL  axis.  The 
resolution  of  this  force  into  A-frame  components  can  be 
accomplished  by  referring  to  Figure  2  where  Euler  angles 
(T ,0 , $)  are  used  to  relate  the  L-frame  and  A-frame  coordinate 
frames : 

Xg  =  -  mg  sin  0  (28) 

=  mg  cos  0  sin  $  (29) 

Zg  =  mg  cos  0  cos  $  (30) 

No  moments  are  produced  by  gravity  because  the  A-frame  origin 
is  located  at  the  aircraft's  center  of  gravity.  Therefore 
Lg,  Mg,  and  Ng  are  zero. 

Non-Gravitational  Forces  and  Moments  -  The  remaining  forces  and 
moments  are  due  primarily  to  aerodynamic,  propulsive  and  control 
effects.  They  are  denoted  simply  X,  Y,  Z,  L,  M,  and  N.  Thus, 
for  example,  Fx  =  Xg  +  X.  These  forces  and  moments  are  developed 
in  Section  V. 

The  A-frame  force  and  moment  equations  (22  through  27)  can 
be  restated  on  the  basis  of  the  above  discussion.  Also  stated 
are  the  relationships  between  A-frame  angular  rates  (P,Q,R)  and 
Euler  angle  rates.  These  relationships  have  been  obtained  from 
Figure  2  by  projecting  the  Euler  rates  (4',0,$)  onto  the  A-frame 
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axes : 


m  [U  +  QW  -  RV  +  g  sin  0]  =  X  (31) 

m  [V  +  RU  -  PW  -  g  cos  0sin$]  =  Y  (32) 

m  [W  +  PV  -  QU  -  g  cos  Geos  4>]=  Z  (33) 

Ix  p  +  (Iz  -  Iy)  QR  -  Jzx  (R  +  PQ)  =  L  (34) 

Iy  Q  +  (Ix  -  Iz)  RP  -  Jxz  (R2  -  P2)  =  M  (35) 

Iz  R  +  (Iy  -  Ix)  PQ  -  Jxz  (P  -  QR)  =  N  (36) 

P  =  i  -  ¥  sin  0  (37) 

•  • 

Q  =  0  cos  $  +  V  cos  0  sin  $  (38) 

#  • 

R  *  ¥  cos  0  cos  $  -  0  sin  4>  (39) 


These  nine  equations  are  an  almost  exact  description  of  the 
motions  of  an  aircraft  operating  near  the  Earth's  surface. 

The  derivation  to  this  point  has  used  only  four  simplifying 
assumptions,  repeated  here: 

1.  Aircraft  mass  is  constant 

2.  The  Earth  can  be  considered  an  inertial  frame 

3.  The  aircraft  is  a  rigid  body 

4.  The  aircraft  is  symmetrical  about  its 
x  -  z  plane. 

The  equations  can  be  further  developed  along  any  one  of  several 
paths.  In  this  case  they  will  be  manipulated  (in  Section  V) 
into  the  form  generally  used  for  linearized  aircraft  stability 
and  control  studies . 
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5.0  Derivation  of  Linearized  Equations  of  Motion 
A  number  of  simplifying  assumptions  are  required  to  develop 
linearized  equations  of  motion  from  the  general  equations 
discussed  in  Section  4.0.  These  approximations  are  (continuing 
the  numbering  sequence  begun  in  Section  4.0): 

Assumption  5:  The  aircraft  is  assumed  initially 
to  be  in  equilibrium  flight  with 
no  linear  or  angular  accelerations, 
no  angular  rates,  and  no  initial  roll 
angle  or  lateral  velocity. 

Assumption  6:  The  X*  axis  is  fixed  in  the  aircraft  and  is 
parallel  to  the  projection  on  the 
XA  -  ZA  plane  of  the  relative  wind 
vector  auring  equilibrium  flight. 

In  other  words,  stability  axes  will 
be  used. 

Assumption  7:  Small  disturbance  (perturbation) 
theory  will  be  used.  Motions  and 
forces  will  be  referred  to  the 
equilibrium  flight  condition  of 
Assumption  1.  Variables  at  this 
flight  condition  will  be  identified 
by  the  subscript  o.  Change  from 
this  condition  will  be  indicated  by 
the  prefix  A.  Thus: 


U  =  UQ  +  AU 

V  =  VQ  +  AV 

W  =  WQ  +  AW 

P  =  PQ  +  AP 

Q  =  Qq  +  AQ 

R  =  R  +  AR 
o 

¥  =  ¥  +  AH' 

o 

0  =  0Q  +  A0 
$  =  $  +  AS- 
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(40) 


Assumption  8:  Small  angle  approximationswill 
be  made.  For  example, 
sin  A0  =  A0 
cos  A0  =  1 

Assumption  9:  Higher  order  terms  will  be  neglected. 

For  example,  (UD  +  AU) (AQ)  -  UQ  AQ 

By  virtue  of  Assumption  5,  P  ,  Qq,  Rq,  4q,  and  VQ  are  zero. 

By  virtue  of  Assumption  6,  WQ  is  zero.  'VQ  is  a  special  case. 

Since  it  is  the  angle  between  the  XL  and  Xq  axes,  it  can  vary 

from  0  to  2 7T .  In  this  analysis,  we  will  specify  that  heading 

angle  be  referred  to  the  initial  Xc  axis.  Thus  VQ  can  be 

treated  as  if  it  were  zero.  The  perturbation  relations  of 

equation  40  can  therefore  be  revised: 

U  =  UQ  +  AU 

V  =  AV 

W  =  AW 

P  =  AP 

Q  =  AQ 

R  =  AR 

¥  =  AY 

0=0  +  A0 

o 

4  =  A4> 

The  derivative  with  respect  to  time  of  each  of  these  variables 
can  be  easily  determined,  e.g., 

*  d 

U  =  _  (U0  +  AU)  =  AU 
dt  ° 
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Equation  41,  Assumptions  8  and  9,  and  fthe  trigonometric 


identities 


sin  (A  +  B)  =  sin  A  cos  B  +  cos  A  sin  B 
i  cos  (A , +  B)  =  cos  A  cos  B  -  sin  A  sin  B  1 

i 

can  be  applied  to  equations  31  through  39  to  produce 


m  [AU  +  g  sin  0O  +  g  cos  0  A0]  =  X 

m  [AV  +  UQ  AR  -  g  cos  0Q  A$1  =  Y 

m  [AW  -  U  AQ  -  g  cos  0Q  +  g  sin  0Q  AG]  =  Z 

I  AP  -  J  AR  =  L 

X'  ZX  I  I 

I  '  AQ  =  M 

y 

I,  AR  -  J  AP.  =  N 
Z  ZX 


P  =  AP  =  A#  -  A^  sin  0C 
Q=  AQ  =  AG 


R  =  AR  =  Af ’ cos  0 
i  o 


Equations  43  can  be  evaluated  at  the  equilibrium  flight  condition 
where  A  quantities  are  zero: 

1  t  ’  I 

mg  sin  0^  =  X„ 

J  o  o 

0  =  Y 


-mg  cos  WQ 


0^  =  i. 


0  =  Lq 

0  -  M_ 


0  =  N, 


Subtracting  these  equations  from  equations  43  produces  the 
following  perturbation  equations: 
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i 


m  [AU  +  g  cos  0  A0]  =  AX 

o 

m  [AV  +  UQ  AR  -  g  cos  0Q  A4>] 


m 

(AW 

-  UG  AQ  + 

I 

AP 

• 

-  J  AR 

X 

zx 

1 

AQ 

=  M 

y 

• 

i 

AR 

i 

c, 

i 

» 

> 

►d 

z 

AX  = 

X  - 

•  XQ,  etc. 

AY 

AZ 


(45) 


where 

Next  to  be  developed  are  the  right-hand  sides  of  equations  45. 

' 

Thie  conventional  practice  of  expressing  these  terms  as  a  Taylor 

j  I 

series  expansion  is  utilized.  The  development  presented  in 
Chapter  4  of'  Reference  2  is  closely  followed. 

The  Taylor  series  is  of  the  form  (neglecting  higher  order  terms) 


AX  =  XA  AA  +  ^  AB+XC  AC  +  .  .  . 


(46) 


where  AA,  AB,  AC  .  .  .  are  the  variables  which  describe  the 
nlotions  of  the  aircraft  or  which  otherwise  contribute  to  AX 

I  i 

and  where 


r  3X 

XA  -  3A  ' 

I 

evaluated  at  the  equilibrium  flight  condition. 


Several  additional  approximations  and  assumptions'  are 

I  * 

desirable  at  this  point: 

Assumption  10:  Higher-order  terms  in  the  Taylor  series 
expansions  can  be  neglected. 

Assumption  11;  The  air  mass  is  fixed  with  respect  to 
the  Earth;  that  is,  there  is  no  wind. 

This  simplifying  assumption  allows 
use  of  components  of  ground  speed  VR 
■ (AU,AV,AW)  as  variables  in  the  Taylor 
series  instead  of  components  of  airspeed. 
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Assumption  12:  Series  terms  involving  the  derivatives 
with  respect  to  these  variables  will  be 
considered:  AU,  AV,  AW,  AW,  AP,  AQ, 

AR,  A6e,  A<5a,  and  A6_.  The  last 
three  variables  are  the  perturbation 
deflections  of  the  elevator,  aileron, 
and  rudder,  respectively. 

Assumption  13:  Series  terms  involving  derivatives  of 

AX,  A2 ,  and  AM  with  respect  to  AV,  AP, 

AR,  A6^  and  A6  can  be  neglected. 

Terms  involving  derivatives  of  AY, 

AL,  and  AN  with  respect  to  AU,  AW,  AW, 

AQ,  anrf  A6_  can  also  be  neglected. 

© 

Assumption  14:  The  derivatives  3X/3A  Q  ,  3X/3A  W  ,AX/A6e  , 
and  AY/A<Sa  are  negligibly  small. 

These  assumptions  are  discussed  and  justified  in  Chapter  4 
of  Reference  2.  Assumption  1!)  is  necessary  in  order  to 
retain  linearity  in  equations  45.  Assumption  13,  while  not 
essential,  is  made  because  experience  has  shown  that  it  is 
a  reasonable  one  (for  most  applications)  and  because  it 
allows  separation  of  the  six  equations  (45)  into  two  sets  of 
three . 

The  definition  of  equation  46  and  the  above  assumptions 
yield  expressions  for  the  non-gravity  perturbation  forces  of 


equations  45: 

AX 

= 

xu 

AU 

+ 

Xw 

AW  + 

AT 

AY 

= 

yv 

AV 

+ 

^AP 

+  Yr 

AR 

+  Y6. 

A6r 

• 

r 

AZ 

= 

Zu 

AU 

+ 

Zw 

AW  + 

Z  . 

w 

AW  + 

Zq  AQ  + 

z«e  «. 

AL 

= 

Lv 

AV 

+ 

Lp 

AP  + 

Lr 

AR  + 

L{  A6a 

cl 

+  L6r  Mr 

AM 

= 

Mu 

AU 

+ 

*w 

AW  + 

*W 

AW  + 

Mg  AQ  + 

M6e 

AN 

= 

Nv 

AV 

+ 

NP 

AP  + 

Nr 

AR  + 

Nr  A  6 

<5a  a 

+  Njr  A6r 

(47) 
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Before  these  expressions  are  substituted  into  equations  45,  one 
question  must  be  resolved:  Equations  45  contain  as  variables 
A-frame  rates  (Ap,AQ,AR)  as  well  as  Euler  angle  displacements 
(At,A0,A4).  A  decision  must  be  made  as  to  which  set  of  variables 
will  be  carried  henceforward. 

For  the  X-Z-M  set  of  equations  45,  there  is  no  difficulty. 

By  virtue  of  equations  44,  Aq  =  AG,  from  which  can  be  obtained 
•  • 

AQ  =  A0.  Thus  AQ  and  AQ  can  be  readily  eliminated  in  favor 
of  A0  and  its  derivatives. 

For  the  Y-L-N  set,  there  is  no  difficulty  either,  except 
that  some  additional  terms  must  be  accepted  for  either  set 
of  variables.  One  can  solve  equations  44  for  AJ>: 

• 

A<t  =  AP  +  tan  0O  AR, 

but  attempts  to  integrate  this  equation  in  order  to  obtain  an 
expression  for  A<1>  (needed  in  the  Y  equation)  encounter  the 
difficulty  that  /AP  and  /AR  cannot  be  uniquely  determined  since 
they  depend  on  the  order  that  rotations  are  taken  about  the 
A-frame  axes.  To  avoid  this  difficulty  (and  to  retain  uniformity 
since  the  Euler  angle  AG  is  used  in  the  X-Z-M  equations)  the 
Euler  angles  A<J>  and  At  are  used  as  the  variables  in  this  analysis. 
Thus  equations  44  are  used  to  eliminate  AP  and  AR  in  equations  45. 

We  will  combine  several  steps  in  arriving  at  the  final 
form  of  the  aircraft  equations  of  motion: 


1.  The  six  equations  will  be  separated  into 
the  X-Z-M  set  and  the  Y-L-N  set. 

2.  AQ  will  be  eliminated  in  fayor  of  A 0, 
and  AP  and  AR  in  favor  of  A$  and  A*!1  by 
using  equations  44. 

• 

3.  The  Laplace  notation  ()  =  s()  will  be 
introduced . 

4.  Equations  (47)  will  be  substituted 
into  equations  (45)  . 

The  X-Z-M  set  -  the  longitudinal  equations  -  are  thus 
developed  from  equations  45  as: 

[ms  -  X  1  AU  +  [-XWJ  AW  +[mg  cos  0O]  A0  =  AT 

[-zu]  AU  +  [ (m-Z^) S  -  Zw]  AW  +  [-(mUQ  +  Zq) s  +  mg  sin  0Q] 

"  [Z6  3  A6e 
2  e 

[-Mu]  AU  +  [-M^S  -  Mw]  AW  +  [IyS  -  Mgs]  A0 

-  tM6e)  4 6e 

TCie  Y-L-N  set  -  the  lateral  equations  -  become: 

[ms  -  Yy]  AV  +  [-YpS  -  mg  cos  0Q]  A4> 

+  [  (mU0  -  Yr)  cos  0O  +  Yp  sin  0Q]  s  AT  =  [Y^^JAfij. 

[-Lv]  AV  +  [Ixs2  -  Lps]  A$ 

+  [-(Ix  sin  0O  +  Jzx  cos  0O) s2  -  (Lr  cos  0O  -  Lp  sin 
=  [LfiaJ  A6a  +  [L6r]  A6r 


(48) 
A0 

(49) 

(50) 


(51) 

0o)s]  AV 

(52) 
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t-Nv]  AV  +  t-Jxz  s2  -  NpS]  A* 

+  [(Iz  cos  0O  +  Jxz  Sin  0o)s2  -  (Nr  cos  0O  -  Np  sin  0O) s]  AY 

=  [Nr  ]  A6  +  [Nr  ]  A6r  (53) 

°a  r 

Equations  48  through  53  are,  to  summarize,  the  six  linearized 
rigid  body  equations  of  motion,  written  in  aircraft  stability 
axes.  Dependent  variables  are  the  perturbed  aircraft  body 
(A-frame)  velocities  (AU,AV,AW)  and  the  perturbed  Euler  angles 
(A0,A4>,AY)  and  their  derivatives.  Longitudinal  (X-Z-M)  and  lateral 
(Y-L-N)  equations  are  not  coupled. 

Coefficients  of  the  variables  in  these  equations  are 
constants  whose  values  are  determined  by  the  aircraft's 
geometric  and  mass  properties  and  its  equilibrium  speed,  UQ. 

These  coefficients  will  be  developed  in  literal  form  in  Section  6.0. 
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6.0  Summary  of  Analytical  Expressions  for  Stability  Derivatives 

The  equations  of  motion  developed  in  Section  V  have  sDecified 
a  set  of  aircraft  stability  derivatives  (e.g.  Xy)  which  are  defined 
analytically  in  this  section. 

Each  derivative  represents  the  (partial)  change  in  a  force 

(or  moment)  due  to  an  incremental  change  in  a  variable  from 

the  reference  value.  Thus  the  derivative  represents  the 

change  in  force  along  the  aircraft  X^  axis  due  to  a  change  in 
forward  speed,  or 


where  the  subscript  o  indicates  that  5X/9U  is  to  be  evaluated 
at  the  reference  (equilibrium)  flight  condition. 

Definition  0f  these  derivatives  here  will  utilize  extensively 
the  material  contained  in  Chapter  4  of  Reference  2.  Emphasis  in 
this  report  will  be  placed  on  assembling  the  desired  material 
into  a  compact  and  usable  form  rather  than  on  repeating  the 
derivations  given  in  the  reference. 

Reference  2  develops  the  derivatives  in  their  non-dimensional 
form.  Table  4.1  of  the  reference  defines  the  relation  between 
these  non-dimensional  coefficients  and  the  dimensional  coefficients 
used  here.  For  example, 


X  =  as  cx  p* 

u  UQ  xu  rps 

Analytical  expressions  for  the  non-dimensional  coefficients  are 

derived  in  References  1,  2,  and  3  and  are  reproduced  here.  For 

example,  on  pp  148-150  of  Reference  2  is  developed  an  expression 

f°r  Cx  : 
u 


Cx  =  -  3  CD  -  C  tan  0Q  (constant  speed  prop) 
u  o  Lo 

Table  I  presents  similar  information  for  each  of  the  required 
derivatives:  Each  derivative  is  defined  in  terms  of  the 

corresponding  non-dimensional  derivative.  Then,  the  expression 
for  the  non-dimensional  derivative  is  given  together  with  the 
source  of  the  expression. 

In  Section  7.0,  these  expression  are  evaluated  for  the 
"Buffalo"  and  "Twin  Otter"  aircraft  at  several  flight  conditions 
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TABLE  I 


ANALYTICAL  EXPRESSIONS  FOR  STABILITY  DERIVATIVES 


V  _  qS  r  lbs 

»  '  U0  xu  fp* 

C  =  -3  C_  -  C  tan  0Q  /constant  J 

xu  o  Lo  °  (speed  prop/ 

(Ref.  2,  pp  148-150) 

x  =  as  c 

W  u0  xa  fps 

Cx  -  CL0  I1-  ireAp)  <Ref-  2-  P  147> 
a  u 

Z  =  qs  c 

u  uc 

C  =  -2  C,  (Ref.  2,  pp  129-130, p  150) 
zu  Lo 

2  =  c  lbs  . 

w  UQ  *a  fps 

c2a  *  -(CLc  +  So*  <Ref.  2,  P  147) 

n  —  ^  _  Q  p  lbS 

w  2U2  q  z  *  ft/sec2 

zuo  a 

p.  S  (tail  contribution) 

c2.=  -2  aT  ^tTs  at  (Ref*  2'  p  165) 
a 

Z  =  qS  C_  — — — 

q  2U0  2q  rad/sec 

l  S  (tail  contributionT 

Czq*-2»T  "T  -§  (Ref.  2,  p  154) 

Z.  =  qS  C,  lb® 

e  6  rad 

e 

da-  S 

%  *  aT  5—  s1  nT  ,Re£-  3'  P  250) 

Se  5e 

qSc  ft  lbs 

Mu  ■  Sr  Cmu  T5T- 

C  =0  (Ref.  2,  p  15l) 

mu 

qSc  ft  lbs 

Mw  =  UQ  Cm  fps 

C  =  CL  (h-hn)  (Ref.  2 ,pl47) 
a  a 

H*  '  — §-  qsc  cn.  £t_ibs_ 

2uo  “  fps 

Am 

C  85  C-  —  (tail  only)  (Ref.  2,  p  165) 
in  •  ^  •  c 

a  a  ^ 

C  _  _  ft  lbs 

Mq  2UQ  q  °  mg  rad/sec 

l  (tail  only) 

Cm  =  C_  (Ref.  2,  p  155) 

niq  ^q  c 

M  =  qSc  C  ft  ll3S 

<5  q  ro*  rad 

e  °e 

C_  =  C  ~  (Ref.  3,  p  250) 

m6  z6  c 
e  e 
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TABLE  I  (continued) 


Y  _  qs  c  lbs 

v  uQ  y 6  fps 


I 


F  do 

-  aF  i*  (1- 


b 

aS 

c 

lbs 

p 

2U 

o 

yp 

rad/sec 

b 

rrC 

c 

lbs 

r 

2uo 

q*3 

yr 

rad/sec 

=  qS 

cv 

lbs 

6 

r 

y«r 

rad 

qSb 

c„ 

ft 

lbs 

'v 

Uo 

c 

p 

fps 

C  -  C  =  -  2  a 
yP  yPF  T 


F 

c  a  Cy  -  2  aF  -g.  — 
r  rF  b 

S  d  olv 

C  =  -  a_  _£  —L 

Ya  f  S  *  d6 


(tail  only) 
(Ref.  2,  p  74) 


(Ref.  3,  p  329] 


%  '  \  +  C‘ 


Bfuselage  3. 


CD  :  See  item  90,  Table  III 


7,,  i  ,  (Ref.  2, 

'SL.  f  =  1.2  /AR  &  •  ^  P  486) 

8  fuse 


C.  =  Cv  ~  (Ref.  2,  p  89) 

*sF  ye  b 


qsb  ct  £LiSi_ 

z,Jo  P  rad/sec 


:  See  item  91,  Table  III 
P 


b  ft  lbs 

-  qSb  Co  - 

2Uq  lx  rad/sec 


cJl  =  CZ  +  c. 

r  r  r 

w  F 


C  =  CT  /4  (Ref.  1,  p  112) 
rw 


C  =  C  Jl  (Ref.  2,  p  175) 


yr  b 


L  =  qSb  C. 
6r  l& 


ft  lbs 


Cp  =  C  .  rr-  (Ref.  3,  p  329) 

*6  y6  b 

r  r 


TABLE  I  (continued) 


=  qSb  C, 


ft  lbs 
rad 


N  =  2Sb 


ft  lbs 


v  UQ  ng  fps 


2aT  *  2 

55—  /  cy  dy  (Ref.  3 ,  p  354) 

y. 


C  =  C„  +  C 

nQ  n  '-n.  .  , 

B  Bp  B  fuselage 

i,p 

c_  =  "  cv  -r-  (Ref.  2,  p  82) 
nBF  % 


n 


B  fuselage 


=  -  1.3 


volume  of  fuseage  h 
Sb  ‘  w 

(Ref.  2,  p  492) 


N  =  —  qSb  C  ft  lbS 


P  2U0  n  rad/sec 


Cn  =  Cn  +  cn 
P  Pw  PF 


n 


n. 


cl  r  1 

«.  _  o  1  1  _  _ a  (App  A,  thisj 

‘pw  4  !  tt  AR_,  report) 

& 

-  C  -1  (Ref.  2,  p  171) 
yPF  b 


b  ^  ft  lbs 

N  =  -  qSb  C 

r  2U 


nr  rad/sec 


Cn  "  Cn  +  cr 


n. 


r 

w 

C 

r  v 


*n  *  -  CD 
"r  w 

w 


(Ref.  1,  p  112) 
(Ref.  1,  p  112) 


Nx  =  qSb  C 


ft  lbs 
rad 


Cn  =  -  C  rf  (Ref.  3,  p  329) 

n  6  ••  6  D 


=  qSb  C 


ft  lbs 
nx  rad 


C  :  no  simple  analytical  expression 


n 


5_  available;  assumed  zero  here. 

a 


7.0  Tabulation  of  Analytically-Determined  Stability  Derivative 
Values  for  DeHavilland  "Buffalo"  and  "Twin  Otter" 

In  this  section,  stability  derivatives  and  other  numerical 
data  required  to  model  these  aircraft  are  presented.  The 
derivatives  are  estimated  using  the  analytical  expressions 
summarized  in  Section  6.0. 

Step-by-step  calculation  of  each  derivative  is  carried 
out  in  Appendix  C.  The  geometric  and  mass  data  needed  in 
Appendix  C  to  calculate  the  derivatives  have  been  obtained 
from  Reference  4;  this  material  is  summarized  in  Figures  3  and 
4  of  this  report. 

For  each  aircraft,  three  flight  conditions  are  investigated. 
These  are:  cruise,  slow  flight,  and  landing  approach.  Table  II 
presents  the  parameters  needed  to  define  each  of  these  flight 
conditions . 

All  stability  derivative  information  is  summarized  in 
Tables  III,  IV,  and  V.  In  Table  III,  the  non-dimensional 
derivatives  calculated  in  Appendix  C  are  collected  and  presented 
for  both  aircraft  at  all  three  flight  conditions.  In  Table  IV, 
the  derivatives  are  presented  in  their  "dimensional"  form, 
utilizing  the  definitions  of  Table  I. 

For  convenience,  the  derivatives  are  also  presented  in 
"normalized"  form  in  Table  V.  The  "normalized"  derivatives  are 
obtained  by  dividing  each  dimensional  derivative  by  the  appropriate 
mass  or  inertia  parameter.  Specifically,  force  derivatives  are 
divided  by  m,  aircraft  mass.  The  roll,  pitch,  and  yaw  moment 

and  I2,  respectively. 


derivatives  are  divided  by  Ix,  Iy, 


I 


TABLE  II 

Definition  of  Flight  Conditions 

i  , 


•  ’ 

‘  i 

Buffalo 

/ 

Twin  Otter 

1  .I, 

Cruise 

BBSs 

Landing 

Approach 

Cruise 

6 

W,  weight,  lbs 

40000 

40000 

40000 

12000 

12000 

Uo,  speed,  fps 

400 

i 

215 

154 

278 

176 

120 

i 

t 

h,  altitude,  ft 

10000 

0 

0 

10000 

0 

0  ' 

Y0,  flight  path 

angle,  .degrees 

'  0 

‘ 

0 

m 

'  '0 

a 

6f,  flap  deflection, 
degrees 

0 

0 

40 

1 

0 

1 

0 

40 

l  i 

1 

i 


Note:  1  Normal  C.G.  used  for  all  flight  conditions. 

2  Max  T.O.  weight  for  "Twin  Otter"  is  listed  as 

11,579  lbs  in  Ref,.  4.  However,  later  editions 
of  Ref.  4  show  an  increase  to  12500  lbs.  The 
weight  to  be  used  here  has  been  arbitrarily 
chosen  as  ^2000  lbs.  ( 

3  0_  =  yQ  in  this  report  because  of  use  of 
stability  axis  system. 
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i  i 

I  i  l 


I 


TABLE  III  Summary  of  Derivatives  in  Non-Dimensional  Form 


Non- 

Buffalo 

Twin  Otter 

Dimensional 

Derivative 

Cruise 

Slow 

Flight 

Landing 

Approach 

Cruise 

Slow 

Flight 

Landing 

Approach 

Cx 

u 

cx 

a 

Cz 

u 

cz 

a 

■  Cz. 
a 

C 

z 

q 

:•••  . 
mu 

Cm 

a ' 

% 

Cm 

q 

c 

mr 

6e 

-.108 

.164 

-.60 

-5.24 

-1.33  ' 

-7.83 

.465 

0 

-.78 

-6.05 

-35.6 

2.12 

-.171 

.441 

-1.54 

-5.26 

-1.33 

-7.83 

i 

.465 

0 

-.78 

-6.05 

-35.6 

2.12 

-.186 

.815 

-2.98 

-5.33 

-1.33 

-7.83 

i  .465 

0 

-.78 

-6.05 

-35.6 

2.12 

-.141 

.234 

-.84 

■  -5.25 

-1.60 

-6.40 

.450 

0 

-.78 

-6.15 

-24.6 

1.73 

-.195 

.435 

-1.56 

-5.27 

-1.60 

-6.40 

.450 

0 

-.78 

-6.15 

-24.6 

1.73 

-.237 

.920 

-3.30 

-5.36 

-1.60 

-6.40 

.450 

0 

-.78 

-6.15 

-24.6 

1.73 

% 

-.362 

-.362 

-.362 

-.492 

-.492 

-.492 

% 

C 

-.055 

-.055 

-.055 

-.085 

-.085 

-.085 

.368 

.368 

.368 

.429 

.429 

.429 

'  Cy<S 

°r 

-.233 

-.233 

-.233 

-.317 

-.317 

-.317 

C* 

6 

-.125 

-.125 

-.125 

-.103 

-.103 

-.103 

c* 

p 

-.53 

-.53 

-.53 

-.60 

-.60 

-.60 

r 

.113 

.231 

.410 

.138 

.233 

.451 

C"6r 

-.024 

-.024 

-.024 

-.024 

-.024 

-.024 

V 

.20 

.20 

.20 

.38 

.38 

.38 

c 

6 

.101 

.101 

.101 

.121 

.121 

.121 

C 

np 

-  .037 

-.134 

-.283 

-.054 

-.129 

-.310 

Cn 

C  r 

-.171 

-.175 

-.188 

-.171 

-.174 

-.191 

n6 

r 

.107 

.107 

.107 

.124 

.124 

.124 

0 

0 

0 

0 

0 

0 
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TABLE  IV  Summary  of  Derivatives  in  Dimensional  Form 


Buffalo  | 

Dimensional 

Derivative 

Cruise 

IS  ^ 

xu 

lbs/fps 

-35.7 

-41.4 

-32.4 

V 

AW 

lbs/fps 

54.3 

106.8 

141.8 

zu 

lbs/fps 

-199 

-372 

-519 

Zw 

lbs/fps 

-1733 

-1273 

-927 

zw 

lbs/fps2 

-5.55 

-7.56 

-7.60 

z 

lbs 

-13,090 

-9560 

-6900 

q 

rad/sec 

2  6 

lbs 

61600 

24150 

12500 

e 

rad 

ft-lbs 

0 

C 

0 

u 

fps 

M 

w 

ft-lbs 

-2610 

-1910 

-1372 

Fps 

\ 

ft-lbs 

-255 

-347 

-347 

fps2 

Ma 

ft-lbs 

6 

- .  6n  Oxio 

-.440xlP6 

-.317xl06 

q 

rad/sec 

«(S 

ft-lbs 

2 . 84xl06 

1J12x106 

,575xl06 

e 

rad 

lbs/fps 

-~.P  Wil 

-87.5 

-63 

Y 

lbs 

-880 

-637 

-460 

P 

rad/sec 

Y 

lbs 

5850 

4260 

3070 

rad/sec 

Y6. 

lbs 

-30,800 

-12,100 

-6,260 

ur 

rad 

Lv 

ft-lbs 

-3980 

-2900 

-2085 

fps 

Lp 

ft-lbs 

6 

-.810x10 

6 

-.591x10 

6 

-.426x10 

rad/sec 

Lr 

ft-lbs 

rad/sec 

+1.7  3xl05 

+  Z58xl05 

+3.30X105 

ft-lbs 

-3 . 06x10 

-1.20x10 

-.62x10 

r 

rad 

6 

6 

6 

ft-lbs 

+2.  55x10 

+1.00x10 

+  .516x10 

a 

_ rad 

Nv 

"“Ft -lbs 
fps 

3215 

2340 

1685 

N 

ft-lbs 

-56500 

-149500 

-22800 

P 

N 

rad/sec 

ft-lbs 

5 

-2.62x10 

5 

-1.95x10 

5 

-1 5 1x10 

r 

rad/sec 

N6 

ft-lbs 

1.365xl06 

.535xl06 

.  276xl06 

r 

rad 

N, 

ft-lbs 

0 

0 

0 

a 

rad 

Twin 

Otter 

- ! 

Cruise 

v  ]\ 

-14.5 

24.1 

-86.5 

-541 

-17.2 

38.2 

-137 

-464 

-14.5 

56.1 

-201 

-327 

-1.93 

-2.16 

-2.19 

-2155 

-1836 

-1253 

12920 

6975 

3265 

0 

0 

0 

-520 

-446 

-306 

-48 

-65 

-65 

-.534x10 5 

5 

-.459x10 

-.3L3xlD5 

.320xL(f 

.17 4x1 06 

.  OSlSxlO6 

-50.6 

-284 

-43.1 

-243 

“  3  0 «  U 
-167 

1430 

1228 

840 

-9100 

-4910 

-2300 

-686 

-589 

-404 

-330000 

-111800 

-76400 

+30000 

+43500 

+57400 

-44500 

-24200 

-11300 

+7.03cl05 

+3.83sd05 

+1.79x1 05 

806 

692 

475 

-11700 

-24100 

-39500 

-37100 

-32400 

-24300 

2.30xl05 

L25xl05 

.585xl<? 

0 

0 

0 
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TABLE  V  Summary  of  Derivatives  in  Normalized  Form 


Normalized 

Derivative 


Twin  Otter 


N  /I 
v  z 

P/;Cz 

r^z 

6x^z 

6a/Iz 


Cruise 


-.160 

-1.397 

-.0045 

-10.53 

49.6 


0 

-.0121 

-.0012 

-2.79 

13.2 


-.097 


-.710 

4.71 

-24.8 


-.0146 

-2.96 

.633 

-1.12 

9.34 


.0072 

-.126 

-.586 

3.05 

0 


-.033 
.086 
-.300 
-1.026 
-.0061 
-7.70 

19.45 

0  “ 

-.0089 

-.0016 

-2.04 

5.17 


-.0705 

-.514 

3.44 

-9.76 

-.0106 

-2.16 

.945 

-.44 

3.66 

.0052 

-.334 

-.436 

1.195 

0 


-.418 

-0.747 

-.0061 

-5.55 

10.08 

0 

-.0064 

-.0016 

-1.47 

2.675 


-.0508 

-.371 

2.48 

-5.05 

-.0076 

-1.56 

1.208 

-.227 

1.89 
TOO 3 8 ~ 

-.510 

-.338 

.617 


Cruise 

Slow 

Flight 

WE$B35Sm 

-.039 

-.0462 

-.039 

.065 

.1027 

.151 

-.232 

-.368 

-.540 

-1.454 

-1.247 

-0.880 

-.0052 

-.0058 

-.0059 

-5.80 

-4.93 

-3.37 

34.8 

18.72 

8.79 

0 

0 

0 

-.0236 

-.0202 

-.0139 

-.0022 

-.0030 

-.0030 

-2.42 

-2.08 

-1.42 

14.52 

7.9 

3.7 

-.136 

-.116 

-.0806 

-.764 

-.654 

-.449 

3.84 

3.30 

2.26 

-24.4 

-13.2 

-6.19 

-.0282 

-.0242 

-.0166 

-5.35 

-4.60 

-3.14 

1.233 

1.790 

2.360 

-1.83 

-.996 

-.465 

29.0 

15.77 

7.36 

.0197 

.0169 

.0116 

-.286 

-.588 

-.965 

-.905 

-.790 

-.593 

5.61 

3.05 

1.43 

0 

0 

0 

It  INERTIAL  COORDINATE  FRAME 

E:  EARTH -CENTERED  COORDINATE  FRAME 

L:  EARTH  LOCAL-VERTICAL  COORDINATE  FRAME 


FIGURE  1:  REFERENCE  COORDINATE  FRAMES 
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L:  EARTH  LOCAL  VERTICAL  COORDINATE  FRAME 
C:  EARTH-AIRCRAFT  CONTROL  COORDINATE  FRAME 
A:  AIRCRAFT  BODY  COORDINATE  FRAME 
EULER  ANGLES 

¥  =  ROTATION  ABOUT  ZL  AXIS 
0  =  ROTATION  ABOUT  Yc  AXIS 
9  =  ROTATION  ABOUT  XA  AXIS 


FIGURE  2:  ADDITIONAL  REFERENCE  COORDINATE  FRAMES 
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DNC4  Mato  hHi  trlipnp  ITOL  aOItty  trtniperl 


DHC-I  BUFFALO 

IHftVrrncni  between  tho  US  and  Canadian 
vrniiHM  arc  a s  follows: 

CV-7A.  l?8  model,  with  2.K30  «**hp  General 
I .!••••« ri«*  TOldK  IO  turl«npni|ki.  Overall  length 
•  7  ft  4  *»i  |23.*i7  m).  DnMigtialmn  may  lie  chang'd 
f.-ll<>uiiii*  transfer  of  rv«<poti«il>ility  for  aircraft 
u.  ll.M  CHt.yory  from  US  Army  to  USAF. 

60>t!L  Canadian  Defence  Force  model,  with 
3,055  eahp  Genera!  Electric  T 04/1*1  turboprom. 
Overall  length  79  A  0  in  (24  08  m).  Otherwise 
similar  to  CV-7A,  with  only  amall  difTercncea  in 
performance. 

WiKOft:  Cantilevar  hiph-wing  monoplane.  Wing 
ecetion  NACA  64»A417-5  (mod)  at  root, 
NACA  83»A815  (mod)  at  tip.  Aspect  ratio 
9-75.  Chord  11  A  91  in  (3  59  m)  at  root. 
5  A  1 1  in  (119  m)  at  tip.  Dihedral  0*  inboard 
of  nunllff.  5*  outboard.  Incidence  2*  30*. 
8wecpl>ark  at  quarter-chord  K  40*.  Con¬ 
ventional  fail-safe  multi-spar  structure  of  high- 
strength  aluminium  allova.  Kull-apan  double  - 
slotted  aluminium  alloy  Aapa,  outboard  sections 
functioning  aa  ailerons.  Aluminium  alloy 
alot-lip  spoilers,  forward  of  inboard  flare,  are 
actuated  by  Jarry  Hydraulic*  unit.  Spoiler* 
eoupled  to  manually -operated  ailerons  for 
lateral  control,  uncoupled  for  symmetrical 
ground  operation.  Electrically -actuated  trim- 
tab  in  starboard  aileron.  Geared  tab  in  oach 
aileron.  Rudder- aileron  interconnect  tab  on 
port  aileron.  Outer  wing  leading-edges  fitted 
with  electrical ly-control led  flush  pneumatic 
rubber  de-icer  boots. 

Fuskuios:  Fail-safe  structure  of  high-strength 
aluminium  alloy.  Cargo  floor  supported  by 
longitudinal  Icoet  members. 

Tail  Unit:  Cantilever  structure  of  high -strength 
aluminium  alloy,  with  fi*ed-incid<  nro  tail  plane 
mounted  al  tip  of  lln.  Elevator  aorodynnmic- 
ally  and  mo** -balanced.  Fore  and  trailing 
serially -hingo< I  rudders  are  |?ourred  by  tandem 
jacks  o|x*rate<|  by  two  independent  hydraulic 
systems  manufactured  by  Jarry  Hydraulics. 
Trim -tab  on  port  elovotor,  spring  tab  on 
starboard  elevator.  Electrically -controlled 

flush  pneumatic  rubbar  de-icer  boot  on  tail- 
plane  loading -edge. 


Landing  Gear:  Retractable  tricyeio  type. 
Hydraulic  retraction,  huhc  unit  aft,  main  linita 
forward.  Jarry  Hydraulics  oleo-pncuiuatio 
shock -alieorbers.  Goodrich  main  wheel*  and 
tyres,  sixe  37  00  x  15*00*12,  pressure  45  lh/sq  m 
(3*10  kgfcuf).  Goodrich  nneo  wheels  and  tvree 
size  8  90  x  12*50,  priwMiiru  38  Ib/sq  in  (2-07 
kg/nnt).  Goodrich  multi-disc  brakes. 


Power  Plant:  Two  General  Electric  T04  turl>o« 
prop  engines  (details  under  entries  for  in¬ 
dividual  version*,  above),  each  driving  a 
Humiltnn  Standard  63E60-I3  thrco-Mudo  pro- 
jwller,  diameter  14  A  0  in  (4*42  in).  Enel  in  on* 
integral  tank  in  each  inner  w  ing,  cajatcity  533 
Imp  gallons  (2,423  litres)  and  ruhlwr  bag  tanks 
in  each  outer  wing,  cu|*u-ity  330  Imp  gallons 
(1,527  litres).  Total  fuel  capacity  1,738  Imp 
gallons  (7.IHHI  litres).  Itufuelling  points  above 
wings  and  in  side  of  fuselage  for  pressure 
refuelling.  Total  oil  capacity  10  Imp  gallons 
(45*6  litres). 


Diversions,  external : 
Wing  span 
Length  oscrr.Ht 
CV-7A 
CC-115 

Height  overall 
Tai  [plane  span 
Wheel  track 
Wheelbase 

Cabin  doors  (each  side): 
Height 
Width 

Height  to  sill 
Emergency  exits  (each 
leading-edge) : 

Height 

Width 

Height  to  sill  approx 
Rear  esrgo  loading  door 
Height 
Width 

Height  to  ramp  hinge 


00  ft  0  in  (20*20  m) 

77  ft  4  in  (23*57  m) 
70  ft  0  in  (24*08  in) 
28  ft  8  in  (8*73  in) 
32  ft  0  in  (9*75  m) 
30  ft  0  in  (0*20  in) 
27  ft  1 1  in  (8*50  m) 

5  ft  0  in  (1*68  m) 
2  ft  0  in  (0*84  m) 
3  ft  10  in  (1*17  m) 
side,  below  wing 

3ft  4  in  (102  m) 

2  ft  2  in  (0*00  m) 

3  ft  Oin  (1*52  m) 
and  romp: 

20  ft  0  in  (0*33  m) 
7  ft  8  in  (2*33  m) 
3ft  10  in  (1*17  m) 


Dimensions,  internal: 

Cabin,  excluding  flight 
length ,  cargo  floor 
Max  width 
Max  height 
Floor  area 
Volume 


deck: 

31  ft  5  in  (9*58  m) 
8  ft  9  in  (2-07  m) 
0  ft  10  in  (2*08  m) 
243  5  »q  ft  (22*83  itf) 
1,713  cu  ft  (48  50  m») 


Areas: 

Wings,  gross  045  sq  ft  (87-8  nf) 

Ailerons  (total)  39  *q  ft  (3*02  irf) 

Trailing-edge  flaps  (total,  including  ailerons) 

280  sq  ft  (20  01  in*) 
Spoilers  (total)  25*2  sq  ft  (2-34  in*) 

Kin  02  sq  ft  (8*55  rrf) 

Rudder,  including  tab  00  sq  ft  (3-57  rrf) 

Tsilplnnc  131*3  sq  ft  (14-07  nfj 

Elevators,  including  tab  81*5  sq  ft  (7*57  n^j 


Wkioiita  and  Loadings: 

0)ierstiiig  weight  empty,  including  3  crew  at 
200  lb  (91  kr)  each,  phis  trapped  fuel  and  oil 
and  full  cargo  handling  equipment 

23,157  lb  (10,505  kg) 


Max  payload 
Max  T-0  weight 
Max  zero -fuel  weight 
.Max  landing  weight 
Max  wing  loading 
Max  powor  loading 


13,843  lb  (6,279  kf) 
41,000  lb  (18.508  kg) 
37,000  lb  (10,783  ky) 
39,000  lb  (17,090  kg) 
43-4  lh/sq  ft  (212  kf/m*) 
7*2  Ib/rshp  (3*27  kg/eshp) 


Performance  (CV-7A,  at  max  T-O  weight): 
Max  level  speed  at  10,000  ft  (3,030  m) 

271  mph  (435  kmh) 
Max  permissible  diving  speed 

334  mph  (537  kmh) 

Max  cruising  speed  at  10,000  ft  (3,050  m) 

271  mph  (435  kmh) 

Econ  cruising  speed  at  10,000  ft  (3,050  m) 

208  mph  (335  kmh) 

*  Stalling  s)>ccd,  *10°  flaps  at  39.000  lb  (17,090  kg) 

AUW  75  mph  (120  kmh) 

:  Stalling  speed,  flaps  up  st  max  AUW 

105  mph  ( 109  kmh) 
Kale  uf  climb  at  S/L  1,890  ft  (575  m)  min 
,  Service  ceiling  30,000  ft  (9,150  m) 

•  Service  ceiling,  one  engine  out 

14,300  ft  (4,300  m) 
T-O  run  on  firm  dry  sod  1,040  ft  (317  m) 

T-0  to  30  ft  (15  m)  from  firm  dry  sod 

1,5*0  ft  (470  m) 

Landing  from  50  ft  (13  m)  on  firm  dry  sod 

1,120  ft  (342  m) 

Lauding  run  on  firm  dry  sod  610  ft  (186  m) 


Source:  Reference  4 


Figure  3 
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DHC-I  TWIN  OTTEft 

Announced  in  Auguil  1084,  the  Twin  Otter  is 
a  STOU  transport  puwi'ml  bv  two  Pratt  «k 
Whitney  (L’Al1)  P'lUA-20  turlmprop  (Migittm. 
I)c*igu  work  began  in  January  liM4.  Construe- 
tinn  of  an  initial  batch  of  Twin  Otters  was  started 
in  November  of  the  an  mo  year  and  the  first  of 
these  flow  ott May  20,  1005. 

At  lie-  beginning  of  1007,  it  total  of  52  Twin 
(liters  had  been  delivered  or  wero  on  order,  with 
options  on  1 1  more.  They  included  eight  for  the 
Clulc.ni  Air  Force,  two  for  Trans- Australian 
An  lines,  mu*  for  the  Canadian  Department  of 
Hands  ami  Forests,  four  for  Aeralpi  of  Italy, 
one  for  Northern  Consolidated  Airlines,  and 
other*  for  Pilgrim  Airlines  ami  Air  Wisconsin, 
l‘SA.  Production  was  scheduled  to  Ik*  at  the 
calf  of  six  a  month  through  1007. 

Cmler  development  for  delivery  in  1008  is  * 
\ it -ion  of  the  Twin  otter  with  more  powerful 
(din  eotipl  Pratt  A  \Mntiiey  I'TtlA  27  turboprop 
engines,  longer  nose  to  provide  more  ha gg»g*» 
sp.ieo,  and  At*W  of  1 2.500  lb  (5,1170  hjt).  The 
following  data  refer  to  the  current  production 

model. 

Tvi*k:  Twin-turboprop  STD  I,  transport. 

Wivt.s;  braced  high-wing  monoplane,  with  a 
-ingle  Ntreuinline-seetion  bracing  strut  on  each 
*nlc.  Wing  section  NACA  0A  series  menu  line; 
NACA  tiolO  (muddied)  thickness  ilistrilmtinn. 
Aspect  lutio  10.  Constant  chord  of  0  fl  6  in 
( I  PS  hi).  Dihedral  2\  Incidence  2J  30'. 
\o  sweephaek.  All-metal  safe-life  MrtirtllK. 
All-metal  ailerons  which  uUo  droop  for  us©  at 
Haps.  Double-dotted  all-metal  full-span 
trading-edge  llaps.  No  spoilers.  Trim  tabs  in 
aileron*.  Pneumatic- bout  de-icing  ci|uipmeitt 
optiunal. 

j’lMJAuK:  Conventional  all-metal  seiui-niono- 
eo<piv  safe-life  atructurc. 

Tut.  Cmt:  Cantilever  all-metal  structure  of 
high  *tivngth  aluminium  alloys.  Fin  integral 
*•  jfh  tutelage.  Fixed-incidence  tedjdane.  Trim- 
tab*  in  rudder  and  port  elevator,  latter  inter¬ 
connected  with  Haps.  Pneumatic  de-icing 
hoots  on  tailplaue  leading-edge  optionul. 

I.vvniso  Gr.iii:  Non* retractable  tricvelo  typo, 
with  steerable  nose. wheel,  Ituhber  shock- 
absorption  on  tiiuiu  units.  ()l**o -pneumatic 
iio-c- wheel  shock-absorber.  l»oodyonr  main 
wh  el  tyres  si/.e  I  M>t>  *  12.  pressure  32  U>V(  ill 
12  25  l»jr  enr).  (Joodyenr  imse-whecl  tyre  si/e 
a  '.m  1 2-5(1,  ore*** tire  31  lb/*cj  in  (2-18  kj/citr), 

Ciiodrieh  hydraulic  brakes.  Provision  for 
alternative  flout  uml  ski  gear. 

pi iw k it  |*i. \\t:  Two  5711  oslip  Pratt  »V  Whitney 
(CAC)  PTdA-20  turboprop  engines,  each 

driving  n  Huit/ell  tin . blade  rev*-rs|h|e.pitch 

lully-fcnt tiering  metal  propeller,  diameter  8  ft 
0  in  (2  1 1  m).  Fuel  m  two  tanks  (8  .  ells)  under 
cabin  Hour;  total  capacity  hi  It  Imp  gallons 
(1.178  litres).  Two  refuelling  points  on  pert 
side  of  fiiM'higc.  (id  capacity  2  Imp  gallons 
(II  litres)  jut  engine.  Fleet ric  de-icing  system 
for  pro|s>||ers  and  air  intakes  optional. 

Accommodation :  Two  scat*  dd* -bv-sidc  <m  flight 
ih-eK.  Neats  for  13- 1 8  p.»****n,:-,rs  iii  inrun 
cabin  Cabin  divi'hd  by  bulkhead  into  naiiti 
pii**enger  or  In  ight  i  oinpai  t  in*  nt  and  baggage 
nr  toilet  < ofipiii tineid.  Ilooi  on  each  side  of 
main  cabin,  nt  i**nr.  lingg.ige  eoiiipariinents  in 
in».*e  i.i.d  aft  rf  cabin,  each  with  ii)  ward- 
hinged  door  on  port  side. 

S\su:\i<;  Hydraulic  *\*t«-in,  pre— lire  I, Ann  lb 
h*»  Hi  (H‘5  Ivl*  em*t.  fnl  Hap*.  brake*  and  no*e. 
wh**e|  M eent::.*.  N«*  pti*  unuitic  *y*t',m.  One 

2 00 A  stnrter-gen-iHlor  on  ouch  •  ngiia  . 

I7».i:i  i it-»vn >  A\n  la/ 1  m*\ii  \t:  Kndio  and  tadar 
to  eu*tni'iel,,-»  ^peeiliention.  lilind  -lit  mg  in*iru- 
mentation  Mandat'd. 

DlMl.V*Ki\s,  |;\UHN\|.i 

Wing  span  t*5  ft  0  tn  ( IH  si  m) 

Length  overall  -111  H  0  in  (15  oil  m) 

H*’’::ht  overall  Is  It  7  in  (.ViiiJ  m) 

Ta i! plane  i> pun  21  ft  (»i:i  ( <»  |H|ii| 

\\  b-  el  track  12  ft  5  m  (3-7h  mi) 

W  heelbase  14  ft  b  ill  (  I  ,Vl  Ml) 

Pn*-eiiger  door  (pott  sid**): 

Height  4  ft  2  in  (I  27  m) 

Witlth  2  ft  0  in  (U-7H  i.i) 

Height  to  *j)l  3ft  IU  in  (1*17  m) 

Pns-rnger  door  (*terb*'uid  *ide) 

H<  ight  3  ft  HJ  in  ( 115  in) 


dc  Havilland  Canada  DHC-6  Twin  0(t«r  twin-turboprop  transport 


.  .  -  •*  MI  Ml  Mi; 

H.-ight  to  Kin  3 f,  inm  (Pi;  ,M) 

Hag:  -age  roniparltiient  door  (im.m  j; 

"••'“I'*  <-  -'ll  -1  II  l»  si  <1  17  in) 

. . .  . nun,  .l.»»  ,■  jm.m  , 

4  d  -  fi  II  l'7  in) 
4  I)  Hill  (M2  |„) 
-ill  :l  fi  Him  (117  m) 


I  MMK\*IO\*.  jvi  KMSAl.: 

Cabin,  excluding  (light  dnl 
..r  toih*t  eonipartnieiit : 
I.emMh 
.»!.«%  *>..*♦(• 

Max  he  e;  I  it 
Floor  at  i  a 
Volume 

Haggage  eompartineiit  (nn- 

Haggage  compartment  (n-m 


k.  galley  Bml  baggage 


Ark.%«: 

Wings,  gro-s 
Aileron-*  (total) 

T railing-edge  flaps  (total) 
Fin 

JSudder,  ilielu-liiig  tab 


lx  ft  t»  in  (5  Ht  m) 
5  ft  3  m  ( I -Mt  in) 
If*  \  1  hi  { I  5(1  in  , 
2  sij  ft  (7-45  m  i 
3h I  nt  ft  (IU  >7  in  ) 

'»*•*)  V  •liiliie 

22  cii  ft  (0  02  lit1) 
ur)  volume 

52  eu  fl  (147  in1) 


420.*<|  ft  (311-02  nr) 
33  2  ft  (3  0S  ) 
112-2  >«,  ft  (10-42  m  l 
4x  (t  *»|  ft  (4  40  in  ) 
■D  u  sej  |i  (3-10  nr) 


'»*«  tiplnne  1 00  mj  ft  («•  i'n  nr) 

Kh-vntors,  inelinliug  tab  35  sij  ft  (3  25  nr) 

Wi.h.ii  rs : 

H»i*ie  opi  rat  mg  a  ♦  ight ,  inebiding  pilot  ( I  7<Mb-= 
77  kg),  in«fto  (fun  lb  45  kg)  and  toll  ml 

0.1 7<»  lb  (2, sou  kg) 

Max  puv  load  (for  Km  mile  —  I  (in  km  range) 

4.4-10  lb  (2.olu  hc> 
Max  T-O  weight  1 1,5 71*  lb  (5.252  kj;) 

Max  landing  weight  1 1.000  Hi  (4.1'ftO  kg ) 

IV  km  nt  vi \si  i:  (at  max  T-f)  weight): 

Max  erni-ing  *pe«d  ut  10,000  h  (3.n  »o  m) 

I  M  inph  (207  kmh) 

Hcnn  rrilNii.*  -pei'il  at  lO.IM'l  ft  (3.0.50  in) 

1 50  inph  (251  l.mli) 
l.amlti.g  -}■«  .  1  04  5  in|.|i  (  Mil  kmh) 

Hate  of  ehnd ,  at  S  h  ),.'».»»!  f.  C72  •■:)  mm 

Ser\  tee  ceiling  25,5no  ft  ^7,7"o  *oi 

Service  reding,  oje-  engine  out  8.5M0  ft  (2.5'Mi  m) 
T-O  t<»  50  ft  (15  in): 

STOI.  1,120  ft  (31 1  m) 

FAH  Pi  3  1 ,700  II  In  in  II.) 

Handing  from  50  ft  (15  m); 

NTDI,  1,021.  ft  (31  I  in) 

CAH  Pt  3  2,100  ft  iimS  m) 

l» niigo  v. itii  max  fuel,  30  min  tesi-rve 

1*20  miles  ( I.4MI  km) 


Source:  Reference  4  Figure  4 
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Appendix  A:  Derivation  of  Analytical  Expressions  for  Selected 
Stability  Derivatives 

Analytical  expressions  for  several  stability  derivatives 
are  not  included  in  available  references  or  are  not  adequately 
developed.  Consequently,  these  needed  expressions  are  developed 
and  presented  in  this  appendix. 


Derivation  of  C„  and  CT, 

-  “  yD  11 

rF 


P 


F 


Rolling  about  the  X  axis  produces  change  in  angle  of  attack 
at  the  vertical  tail.  This  angle  change  produces  a  side 
force  at  the  tail  and  a  yawing  moment  about  the  axis. 

The  restoring  rolling  moment  produced  is  negligible  compared 
to  the  wing's  contribution. 


From  the  sketch,  the  local  surface  angle  of  attack  is 


0  p  < z)  * 


Side  force  generated  is  then 


dYp  =  -  I  p  UQ2  aF  *  ^  CF(  z)  dz 
^  o 


(assuming  constant  lift  curve  slope 
aF  on  fin  elements) 


dz 


A2 


Assume  that 

cf(z)  -  cr 


Then 


Y 

F 


P 

aF  DJ 


bF 

I 


— 2 - 


c  b  2 
r  F 


6 


1  + 


Next,  define 


3  /c 
t'  r ' 


and  approximate  fin  area  as 

+  ct\  crbp 


SF  “  bP 


-  (1  +  ^p) 


So 


1  _  ~  2  .  P _ l 


V  p  V  a 


F  U  3  F 
o 


1  + 
IT 


and 


fF/*s 


(Pb/2U0) 


?  S  b 

-  -  a  F  - 

3  FS"  b 


Similarly, 

C. 


n 


-  C 


% 

F 

y  "S' 
Yp 
F 


2  bF  SF 

^  7 


3  F  b 


A3 


I 


A  roll  rate  (p)  produces  a  differential  change  in  angle  of 

attack  (a)  on  each  wing.  Thechange  in  a  produces  a  differential 

.  '  .  1 

inclination  and  a  differential  change  in  magnitude  of  the  lift 

and  drag  vectors  on  each  wing,  with  a  resulting  yawing  moment. 


I 


The  local  angle  of  attack  at  a  spanwise  station  y  is 

,  i  .•  i  py 

“(y)  =  wo  +  u 

o 

The  lift  and  drag  on  a  chordwise  element  of  width  dy  are 

f 

i 

1  2 

dL(y)  =  -  p  U.  c(y)  a  n(y)  dy 
2  ° 

and 

12  2 
dD (y)  =  -  p  uQ  c  (y )  [CD  +  k  (aa(y))  ]  dy  k 

^w 

The  yawing  moment  on  a  chordwise  element  of  width  dy  at 
spanwise  station  y  is 

dN(y)  =  -  dL(y)  •  a(y)  •  y  +  dD(y)  •  y 

or 

dN  (y) 

where 

Now  calculate  yawing  moment: 

b/2 

N  =  /  dN  (y) 

-b/2 


=  q  [-a (a.  + 


EZ 


)  +  C 


Df 


+  ka  (a_  + 


EL) 
u  J 


w 


] c (y) ydy 


q  =  -  p  U  2 
H  2  0 


1 

ttAR 


A5 


So 


N 


^2  2  D  v  ^  2  3 

7a  y  +  2ka^  aQ  “-  +  ka  /H-*-  J  dy 

o  /  U0^ 


+  Cr 


y  +  k  a‘ 


w 


Terms  marked  f will,  when  integrated,  equal  zero  by  symmetry, 
so 

b/2 

N  =  [-2  a  otQ  q  +  2  k  a2  a  £-  q]  /  y2  c(y)  dy 
0  o  -b/2 


Next,  evaluate  integral  assuming  an  elliptical  wing  planform 
and  obtain,  for  the  integral, 


*  JL  ,b  3 

T  2  J 


1  fV  q 

4  (2)  S 


So 


N  *  - 


1  ,  .  2  £b 

-  (a  aQ  -  ka  a  ) 


<J3b 


or,  since  CL  =  a  aQ, 
o 


N  - - -  (1-  -~)  •  qSb 

4  ttAR  2Ud  ^ 


Finally , 


N/q  Sb 


nD  (Pb/2UQ) 


o  a 

r  u  -  as1 


w 


A6 


Appendix  B.  Estimate  of  Moments  and  Products  of  Inertia 


No  moment  and  product  of  inertia  data  have  been  located 
for  the  "Twin  Otter".  Reference  5  gives  data  for  the  "Buffalo" 
but  does  not  specify  the  axes  used.  Accordingly,  estimates 
have  been  made  as  necessary.  The  basis  of  these  estimates  is 
presented  in  this  Appendix. 

In  this  report,  the  following  values  of  the  inertia 
parameters  will  be  used: 


Buffalo 

Twin  Otter 

Ix,  sl-ft2 

273000 

24300 

Iy,  sl-ft2 

215000 

22000 

Iz,  sl-ft2 

447000 

41000 

J  ,  sl-ft2 

0 

0 

xz 

f  - - 

i 

Notes : 

! 

1  These  values  will  be  assumed  to  apply 
to  the  A-frame  axis  system.  It  will 
be  assumed  that  changes  between  flight 
conditions  are  negligible. 

2.  Although  Reference  5  gives  a  value 
of  Jxz »  it  does  not  specify  the  axis 

system  used.  This  product  of  inertia 
term  is  very  sensitive  to  small  angular 
changes  in  the  axis  system  orientation. 
It  is  considered  advisable  to  assume 
J  =0  rather  than  use  a  value  whose 
basis  is  unknown. 
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MOMENTS  OF  INERTIA 


1.  DHC- 


Note 


2.  DHC 


"Buffalo"  -  Reference  #5  gives  moments  of  inertia 
for  a  gross  weight  of  38,000  lbs. 


Ix  =  273,000  slug-ft2 

I  =  215,000  slug-ft2 

Iz  =  447,000  slug-ft2 

J  =  +36,000  slug-ft2 

A  Z 


Axis  system  unspecified  in  Reference  5.  Waterline 
axes  assumed. 


Twin  Otter"  -  Corresponding  data  was  not  available 
for  this  aircraft,  and  the  moments  of  inertia  were 
estimated . 

I  =  24,300  slug-ft2 

Iy  -  22,010  slug-ft2 

Iz  =  41,020  slug-ft2 

The  weight  of  each  item  in  the  aircraft  was  evaluated 
using  information  in  Sechler  &  Dunn,  Airplane  Structural 
Analysis  &  Design,  and  Reference  4.  The  mass  of  each 
component  was  then  calculated,  and  used  in  the  appro¬ 
priate  formula  to  determine  moment  of  inertia  in  a 
waterline  axis  system,  as  shown  in  the  following. 
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DHC  -  "TWIN  OTTER" 
ESTIMATED  WEIGHT  BREAKDOWN 


Item 


Weight 
W,  lbs 


Mass  , 
m,  slugs 


Method  of  Estimation 


fuselage 

wing 

ailerons  & 
flaps 

elevators 

horizontal 

stab. 

rudder 

vertical 

fin 

fuel  tanks 

engines 

engine 

installation 

props 

main  gear 

nose  gear 

radio 

instruments 

furnishings 

miscellaneous 

pilots 

oil 


1000 

1500 

50 

35 

130 

45 

155 

150 

600 

200 

350 

250 

150 

100 

100 

500 

575 

340 

40 


1 


31.0 

48.0 

5.13 

6.21 

4.65 

24.8 

10.9 
7.76 

4.66 
3.1 
3.1 

15.5 

17.8 

10.55 

1.24 


8.3%  X  G.W. 
12.5%  x  G.W. 


wing  =  1550# 


•85  lb/ft2  x  56  ft2  , 
1.0  lb/ft2  x  35  ft2 


2.0  lb/ft2  x  65  ft2 
1.3  lb/ft2  x  34  ft2 


horizontal  tail  =  165# 


vertical  tail  =  200# 


2.3  lb/ft2  x  66  ft2 

.48  lb/gal  x  315  imp.  gal. 
lbs 


289 


/engine  -  Jane's 


engines  =  800# 


1.7%  x  G.W. 
3%  x  G.W. 

2%  x  G.W. 
1.25%  x  G.W. 


~  25  lbs/pass. 


2  x  170  lbs 


lb 


empty  wt 

6270 

194.4 

passengers 

3230 

100.2 

fuel 

2500 

77.6 

gross  weight 

12,000 

372.2 

5  gal.  x  8  /gal. 


19  persons  x  170  *b/person 

8  ^b/gal  x  315  imperial  gal. 

method  of  wciqht  estimation  from: 

E.  Scchlcr  &  G.  Dunn. 

Airplane  Structure]  Analysis  &  Design 
Dover  Publications,  New  York,  1963 
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DHC  TWIN  OTTER 


MOMENT  OF  INERTIA  ABOUT  THE  LONGITUDINAL  AXIS 


EQUATION 


VALUES 


INERTIA 


m  +  z2} 

12 

b  =  65.0'  Z  -  -3.5' 

1.745  x  10 

m  {y2  +  Z2} 

y  =  9.3'  Z  =  -3.5' 

.246  x  10 

m  {y2  +  Z2} 

y  =  9.3'  Z  =  -3.5* 

.108  x  10 

m  •  R2 

R  =  3.5' 

.038  x  10 

m  •  R2 

R  =  3.5' 

.019  x  10 

m  •  Z2 

Z  =  5.4' 

.018  x  10 

m  {ki  +  Z2} 

12 

b  =  21.2*  Z  =  -4.5' 

.030  x  10 

m  •  R2 

R  =  4.0' 

.132  x  10 

m  {y2  +  Z2} 

y  -  2.0'  Z  =  -1.0' 

.050  x  10 

m  •  R2 

R  =  7.0' 

.038  x  10' 

m  •  Z2 

Z  =  -4.5' 

.0095  x  10' 

Note : 

components  not 
shown  are  a$umed 
to  contribute 
negligible  Ix 


»« 

4 


>4 

,4 

.4 


>4 

>4 

.4 


2.434  x  10 
2 


I 


slug-ft 


l 


DHC  TWIN  OTTER 


MOMENT  OF  INERTIA  ABOUT  THE  LATERAL  AXIS 


ITEM 

EQUATION 

VALUES 

INERTIA 

fuselage 

\ 

pilots 

n{ (§_  +  *  )  +  x2} 
2  12 

2, 2  j 

R=3 . 5  '  £=41 '  x=-l . 5  ' 

.460  x  104 

4 

passengers , 
furrdAings  j 
nisc. 

n{l2  +  x  } 

„2 

£  =  16. 51  x  =  3.35' 

.480  x  10 

fuel  &  tanks 

£  y 

m{15  +  x  1 

£  =  8.0'  x  =  3.5' 

.144  x  104 

wing 

.2  ;  ! 
mll2  +  <f>  > 

£  =  6.5'  c  =  6.5' 

.029  x  104 

engines 

m  •  R2 

R  =  5.0' 

.062  x  104 

props 

m  •  R2 

R  =  7.0' 

.053  x  104 

main  gear 

m  •  R2 

R  =  4.5' 

.016  x  104 

nose  gear 

m  •  R2 

R  =  13.0' 

.125  x  104 

horiz. tail 

m  .  R2 

R  =  25.8' 

.342  x  104 

vert,  tai] 

m  (x2  +  Z2) 

V 

Note : 

Components  not 
shown  are  assumed 
to  contribute 
negligible  I 

x  =  -27.5'  Z  =  -5.4' 

.490  x  104 

2.201  x  104 

2 

slug-ft 

B5 


. 

■i 

. 

i 


DHC  TWIN  OTTEP. 


MOMENT  OF  INERTIA  ABOUT  THE  VERTICAL  AXIS 


ITEM 


EQUATION 

VALUES 

INERTIA 

I7  z  I 

L  _  v 

Fus  JFus 

.460  x  104 

I  =  I„ 

.480  x  104 

z  y 

2 

l 

4 

m  •  - 

l  =  8.0 

.044  x  10 

m  {(J)2  +  ^  + 

12 

"  2 
<5»  1 

c  -  6.5*  b  =  65' 

1.720  x  104 

m  {y2  +  x2} 

y  =  9.3'  x  =  5.0' 

.278  x  104 

m  {y2  +  x2} 

y  «  9.3*  x  =  7.0' 

.148  x  104 

m  •  y2 

y  =  6.5' 

.033  x  104 

m  •  x2 

x  =  13.0* 

.125  x  104 

m  •  x2 

x  =  -25.8' 

.342  x  104 

2 

m  •  x* 

x  =  -27.5’ 

.472  x  104 

4.102  x  104 

slug-ft2 

Note : 

components  not 

shown  are  assumed 

to  contribute 

negligible  Iz 

fuselage 

pilots 

passengers 

furnishings 

misc. 


wing 

engines 
props 
main  gear 
nose  gear 
horiz.  tail 
vert,  tail 


Derivation  of  Expressions  for  Change  in  Inertias  with 
Rotation  of  Axes 


Inertia  values  are  dependent  on  the  set  of  axes  chosen. 
In  this  part  of  Appendix  B,  equations  are  derived  which 
show  the  effect  on  inertia  values  produced  by  rotating  the 
reference  axes  through  an  angle  £  about  the  Y^  axis.  (The 
Ya  axis  is  universally  chosen  to  be  normal  to  the  aircraft's 
plane  of  symmetry,  so  reorientations  of  the  XA  and  ZA  axes 
only  need  be  considered.) 


From  the  sketch,  the  moments  and  product  of  inertia  in  the 
X  Y  Z  axis  system  are: 

ix  =  /  (y2  +  z2)  dm 

Iy  =  /(x2  +  z2)  dm  (Bl) 

I2  =  /(x2  +  y2)  dm 

J  =  +  /  x  z  dm 
xz  1 
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I 


and,  in  the  X  Y  z  axis  system,  they  are: 


V 


xz 


-  / 

(y-2 

+  z'2) 

dm 

-  J 

(x-2 

+  z-2) 

dm 

-  / 

(x-2 

+  y'2) 

dm 

=  / 

* 

X  Z 

dm 

(B2) 


Also  from  the  sketch,  the  following  relations  can  be 

/  /  ;  *  t 

obtained: 

I 

I 

x '=  x  cos  e  -  z  sin  e 

*  { 

.  y  -  y 

z'  -  z  cor  e  +  x  sin  e 

•'  t  r  t 

.  i  • 

The  desired  relations  are  obtaihed  by  substituting  equations 
(B3)  into  equations  (B2)  and  by  making  use  of  the  trigonometric 
identity 


(B3) 

i 


2  2 
sin  e  +  cos  e  =  1 


These  relations  are 


I  ' 
x 


2  2 

Ix  cos  e  +  Iz  sin  e  +  2  Jxz  cos  e  sin  e 


I. 


I  " 
z 


2  2  ' 

J_  cos  £  +  I  sin  e  -  2  J  „  cos  e  sin  e 

Z  X  X«i 


(B4) 


J  =  J  'cos  2e  +  (I_  -  I  )  cos  e  sin  e 
xz  xz  z  x  I 


I  I 


I  t 
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-  ^  • 


I  I 

-  i^ir— 


i  t 


Estimation  of  Angle  e  Between  X.  .  and  X.  Axes 

J  *  waterlvne  A 


In  this  part  of  Appendix  B,  an  estimate  is  made  of  the 
magnitude  of  the  angle  e  that  is  associated  with  the  use 
of  stability  axes  for  the  three  flight  conditions  being 

I 

considered. 


XW/L  ^waterline 
(parallel  to  fuse 
ref  line) 


e  (negative  as  shown) 


(parallel  to  equilibrium 
remote  wind) 


Assume  that  wing  incidence  has  been  chosen  by  the  manufacturer 
to  produce  a  level  fuselage  attitude  when  the  aircraft  is  in 
flight  at  "Economy  Cruise  Speed"  at  10,000  ft  and  at  an 
arbitrarily-chosen  average  gross  weight, 


For  Buffalo: 


W  =  W  -  M-Payload 

average  T.O  max  2 


=  41000  -  13843/2  (Reference  4  or 

Figure  3,  this  report) 

=  34078  lbs 


VEcon  cr  ~  208  mph  =  306  fps  ®  10,000  ft  (same) 

qEc  =  \  p  UQ2  s,-002754  <306>2  =  82*°  lbs/ft2 


CT  =  W/qS  =  34078/(82)  045)  =  .44 
"EC 
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i 


For  Twin  Otter: 


W 

average 


11/579 


4430 


=  9364  lbs 


(Reference  4  or 
Figure  4,  this  report) 


VEcon  cr  =  156  mPh  =  230  fps  @  10,000  ft  (same) 


.001754  2  2 

•  (230)  =  46.5  Ibs/ft 


CL  =  9364/(46.5) (420)  =  .48 
Ec 


Next,  assuming  slope  of  lift  curve  is  as  calculated  in  Section 
VII  (=  5.2/rad  for  both  aircraft),  one  can  sketch  lift  curves: 


From  these  sketches,  one  can  infer  e  (angle  between  XA  axis 
and  Xjj/l  axis)  at  a  given  CL: 


e 


5.2 


(rad) 
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Now,  calculate  this  angle  for  each  flight  condition  considered 
in  this  report: 


Aircraft 

Flight 
Cond . 

Cl0 

(Section  VII) 

CL 

Ec 

-e 

rad 

deg 

Buffalo 

Cruise 

.30 

.4 

4 

-.0270 

-1.55 

Slow  Fit 

.77 

.0635 

3.64 

Lndg  Appr 

1.49 

.202 

11.6 

Twin  Otter 

Cruise 

.42 

.48 

-.0115 

-.66 

Slow  Fit 

.78 

.0577 

3.30 

Lndg  Appr 

1.65 

.225 

12.9 

Note:  Calculation  has  not  accounted  for  flap 
effects.  Use  of  flaps  (as  in  Landing 
Approach  Flight  Condition)  will  reduce 
|e|  substantially. 


The  above  table  shows  that  |£|  is  largest  for  the  Landing  Approach 
flight  condition.  Accordingly#  the  change  in  inertia  magnitude 
due  to  axis  system  rotation  will  be  investigated  at  this  flight 
condition . 

Equations  (B4)  and  the  waterline-referenced  inertias  of  page  B1 
are  used  to  prepare  the  following  table  (assuming  Jxz  =  0) : 


Buffalo 

Twin  Otter 

(e  =>  -11.6°) 

(e  -  -12.9° 

r*  -  xw swt2 

273000 

24300 

I  '  -  Ix  ,  Sl-ft2 
x  A 

280000 

25040 

%  change 

2.5% 

3% 

I  =  I  ,  sl-ft2 

447000 

41000 

2  ZW/L 

V  -  J»  >  sl-ft2 

440500 

40135 

A 

%  change 

-1.5% 

2, 

This  table  shows  that  inertia  values  are  only  slightly  affected 
by  small  axis  rotations.  It  is  concluded  on  the  basis  of  this 
table  that  variations  in  inertias  due  to  small  axis  rotations 
can  be  reasonably  neglected. 
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Appendix  C  Calculation  of  Stability  Derivatives 


This  appendix  shows  the  step-by-step  calculations  carried 
out  to  obtain  the  derivatives  presented  in  Section  VII. 

In  this  appendix,  the  parameters  needed  to  evaluate  the 
non-dimensional  stability  derivatives  are  developed  (Items  1  -  72) . 
The  source  of  each  parameter  is  shown.  Then,  the  non-dimensional 
stability  derivatives  are  calculated  (Items  73  -  99)  using  the 
expressions  presented  in  Section  VI,  Table  I.  Finally,  in 
Items  100  -  111,  certain  dynamic  constants  (e.  g.,  c/2U0)  needed 
to  calculate  the  dimensional  derivatives  are  tabulated.  Also 
presented  is  a  summary  of  moments  and  products  of  inertia 
developed  in  Appendix  B. 

Some  of  the  parameters  appearing  in  this  appendix  are  used 
only  here,  and  are  not  included  in  the  list  of  symbols  appearing 
in  the  front  of  this  report.  However,  these  parameters  are 
fully  defined  where  they  are  introduced  in  Appendix  C. 


Appendix  C 

Calculation  of  Derivatives 


C2 


C3 


NO. 

PARAMETER 

SOURCE 

VALUE  '  | 

BUFFALO 

TWIN  OTTER 

34 

C  ,  parasite  drage  coefficient 

r  of  aircraft  =  f/S 

(def .) 

.032 

.039 

Note:  These  values  of  CE  were  i 

ised  to  ca! 

.culate  ma) 

.  forward 

speeds  and  rates  uf  of  cl: 

.mb  to  comj 

lare  with  1 

igures 

published  in  Reference  4. 

The  agre< 

iment  was  c 

ood . 

Analytical  estimates  of  C 

(or  f)  s\ 

lould  only 

be 

considered  to  be  rough 

’f  estima1 

:es.  The  \ 

alues 

shown  here  are  reasonable 

ones ,  adec 

[uate  for  1 

his 

study's  purpose. 

35 

p,  atmospheric  air  density, 

Table  11^ 

slugs/ft3 

Cruise 

Std 

.001756 

.001756 

Slow  Fit 

^mos 

.  002378 

.002378 

Lndg  Appr 

Tables 

.002378 

.002378 

36 

U  ,  aircraft  forward  speed 

Table  II 

at  equilibrium  flight 

condition,  fps 

Cruise 

400 

278 

Slow  Fit 

215 

176 

Lndg  Appr 

154 

120 

37 

2 

q,  dynamic  pressure,  lbs/ft 

1  2 
“ipD 

2  o 

Cruise 

140.5 

68.0 

Slow  Fit 

- 

55.0 

36.9 

Lndg  Appr 

28.4 

17.3 

38 

CT  ,  aircraft  lift  coefficient 

o  at  equilibrium  flight 

condition  =  W/qS 

Cruise 

.30 

.42 

Slow  Fit 

- 

.77 

.78 

Lndg  Appr 

1.49 

1.65 

C7 


r  VALUE  | 

BUFFALO 

TWIN  OTTER 

NO 


PARAMETER 


SOURCE 


39 


40 


CD  ,  aircraft  drag  coef.  at 

o  equilibrium  flight  condition 
2 

=  C_  +  CT  /reAR 
Df  '  Lo 
,  Cruise 
Slow  Fit 
Lndg  Appr  , 

I  »  i 

0o,  angle  of  rotation,  at 

equilibrium  flight  ‘  condition, 
of  axis  from  horizontal 

plane,  rad  or;  deg 

Cruise 
Slow  Fit 
Lndg  Appr 


Note:  Item  39:  Landing  Approach 

1  However,  insufficient  inf$ 
1  geometry  to  make  a  valid 
,  the  "noflap"  <?ata  is  appl 
of  the  required  flap  data 


Condition 
rmation  is 
Estimate  o: : 
ed  pending 


.036 

.057 

.127 


.047 

.065 

.155 


Fig.  2, 
this  rpt1. 


0 
0 

-7.5°, 

-.131  red  J-.131  rad 


0 

0 

-7.5°. 


involves  use  of  flaps, 
availably  on  flap  1 
AC  I  Therefore, 

“flaps .{  receipt 


I 


NO .  I  PARAMETER 


41 


-(h.-!^)  ,  stick-fixed  static 

margin,  fraction  of  c 

Note : 

It  will  be  assumed  in  this 
report  that«  aircraft  parameters 
are  such  that  the  static  margir 
is  limited  to  values  between 
.25  (most  fdw  C.G.)  and  .05 
(most  aft  C.G.) 

Therefore: 

(h-h  )  (most  fwd  C.G.) 
n  (normal  C.G.) 

(most  aft  C.G.) 

In  this  report,  the  normal  C.G. 
value  will  be  used. 


_ VALUE _ 

SOURCE  BUFFALO  |  TWIN  OTTER 


Ref.  2, 
Ref.  3 


-.25 

-.15 

-.05 


-.25 

-.15 

-.05 


h 

*  n 

he 

- - 

4 - 

center 

of 

gravity 


stick -fixed 

neutral  point  »  aerodynamic  cente 


(not  used) 


45 


NO. 


PARAMETER 


SOURCE 


VALUE 


BUFFALO 


TW.TN  OTTER 


46 

47 

48 

49 

50 


S  /vertical  tail  area  (including 
F  rudder) ,  ft2 


b  ,  vertical  tail  span,  ft 
F 


,  vertical  tail  chord, ft 


geom 


AR  ,  aspect  ratio  of  vertical 

F  2 
fin  =  bF  /Sp 


AR  ,  effective  aspect  ratio  of 

Fe  vertical  fin  (to  account 

for  stabilizer  end-plate 

effect)  =  1.55  AR_ 

F 


Effective  aspact  ratio  vertical  tall,  At 

Figciib  8-8.  Slope  of  lift  curve,  verti¬ 
cal  tail.  From  NACA  TN  775, 
‘'Analysis  of  Wind-tunnel  Data  or. 
Directional  Stability  and  Control,”  by 

II.  K.  Paaa. 


Fig.  3,4 
Fig.  3,4 
Fig.  3,4 


(def.) 


Ref.  3, 
p  325 


•*r 

CM 

m 

Q, 

ro 


4-1 

<U 

« 


{ 


152 

14.0 

3. 3 (root) 
8.  5 (tip) 


1.29 


2.00 


82 

11.0 

9.2  (root) 
5.7  (tip) 


1.48 


2.30 


CIO 


valu£  I 

NO. 

PARAMETER 

SOURCE 

BUFFALO 

TWIN  OTTER! 

51 

a  ,  lift  curve  slope  of  vertical 
P 

tail,  deg-^  or  rad 

above 

figure 

/  .044/deg 
\ 2 . 5/rad 

.049/deg 

2 . 8/rad 

52 

A_,  taper  ratio  of  vertical  tail 

F  =  c_  /c 

P  '  p 

tip  root 

Item 

.64 

.62 

53 

£  ,  distance,  wing  quarter  chord 
F  to  vertical  tail  quarter 
chord,  ft 

Fig.  3,4 

44 

25.5 

54 

Z_,  vertical  distance,  aircraft 
CG  to  mean  chord  of  vertical 
tail,  ft 

Fig.  3,4 

10 

5 

55 

dd/dfl,  rate  of  change  of  side- 
wash  angle  at  vertical 
tail  with  change  in  side¬ 
slip  angle,  rad/rad. 
Estimate  requires  wind 
tunnel  data.  Will 
arbitrarily  use  value  of 
.1  here. 

(Est) 

.1 

.1 

56 

2 

Sr ,  area  of  rudder,  ft 

Fig.  3,4 

60 

34 

57 

sr/sF 

'54*46 

.4 

.4 

58 

dot  /dfi-.,  rudder  effectiveness 

F  r 

chart  @ 
Item  22 

.58 

.58 

59 

Z  ,  vertical  distance  of  root 
w  quarter  chord  below  fuselage 
center  line,  ft 

Fig.  3,4 

-3.5 

-3.0 

60 

h,  fuselage  height  at  wing  root, 
ft 

Fig.  3,4 

9.8 

7.0 

61 

w,  fuselage  width  at  wing  root, 
ft 

Fig.  3,4 

10.0 

6.0 

Cll 


1 


J. 


PARAMETER 


VALUE 

SOURCE  BUFFALO  ]  TWIN  OTTER 


ca/c ,  ratio,  aileron  chord  to 

local  wing  chord  Fig.  3,4 


man 


lav 


Figure  9-16.  Aileron  effectiveness. 


t,  aileron  effectiveness 


above 

figure 


y„  distance,  aircraft  centerline 

to  inboard  end  of  aileron, ft  Fig.  3,4 

y_,  distance,  aircraft 

centerline  to  outboard  end 
of  aileron,  ft  Fig.  3,4 

3 

Volume  of  fuselage,  ft  Fig.  3,4 


14.0 


31.0 


1300 


i 


'  7. 


NO .  PARAMETER 


SOURCE  | BUFFALO 


TWIN  OTTER 


67  Cq  ,  wing  profile  drag 
'  I  f  coefficient  * 


C  ,  wing  drag  coefficient* 

=  cn  +  CT  2/ttA 
fw 


Cruise 
Slow  Fit 
Lndg  Appr 


Ref.  Cl 
(below) 


Ref.  3, 
P  73 


.006  .006 

(representative  values 
from  data  plots  of 
ref.) 


.009 

.025 

.077 


Reference  Cl:  "Theory  cf  Wing  Sections  Including  a 
Summary  of  Airfoil  Data",  Abbott 
&  Von  Dcenhoff,  Dover  Publications, 
Inc. ,  N, Y. ,  1959. 


*  See  note  at  Item  (39 


(not  used) 


C13 


NO. 


BUFFALO 


yn — @- 


TWIN  OTTER 


73 


°,u  *  -  3CDc  -  \  tan 


Cruise 

=  -  3  ( .036)  -  0  =  -  .108 
Slow  Fit 

=  -  3 ( .057)  -  0  =  -  .171 
Lndg  Appr 

=  -  3 ( . 127)  -  (1.49) (-.131) 
=  -  .186 


Cruise 

=  -  3  ( .047) 
Slow  Fit 

=  -  3  ( . 065) 
Lndg  Appr 
=  -  3 ( . 151) 


-  0  =  -  .141 

-  0  =  -  .195 

-  (1.65) (-.131) 
=  -  .237 


74 


a 


1  - 


Clo  11 


=  .547  C 


fr- 


2  =l£> 

7re  AR 

®  &  , 

2(5.2)  ^ 

7 rt .  ^  5 )  (9.75 )  J 


=  CT 


I"  2(5.2)  I 

L1  "  TT  ( -75)  (10)  J 


Cruise 

=  .547  (.30)  =  .164 
Slow  Fit 

=  .547  (.77)  =  .441 
Lndg  Appr 

=  .547  (1.49)  =  .815 


=  .558  C 

Lo 

Cruise 

=  (.558)  (.42)  -  .234 
Slow  Fit 

=  (.558)  (.78)  =  .435 
Lndg  Appr 

=  (.558)  (1.65)  =.920 


C14 


NO . 

BUFFALO 

TWIN  OTTLR 

75 

38 

C7  "  "  2  CL 

u  Lo 

Cruise 

Cruise 

=  -  2  ( .30)  =  -  .60 

=  -  2  (.42)  =  -  .84 

Slow  Fit 

Slow  Fit 

=  -  2  ( .77)  =  -  1.54 

=  -  2  (.78)  =  -  1.56 

Lndg  Appr 

Lndg  Appr 

=  -  2  (1.49)  =  -  2.98 

=  -  2  (1.65)  =  -  3.30 

76 

Define  V  ,  horizontal 

T 

tail  volume: 

10  9 

ST 

v  =  n  —  — 

T  T  c  S 

14  ^3  N  1 

(46)  (233) 

(25)  (100) 

(1,U)  (10.1)  (945)  "  1,12 

11 '  ;  (6.5)  (420) 

.  915 

77 

C  =  -  2  am  Vm 

z .  -  T  T  da 

11  76  19 

— __J 

=  -  2(3.5)  (1.12)  (  .17)  *  -  1.3C 

=  -  2(3.5)  (  .915)  ( .25)  = 

-  L60 

C15 


i 


NO. 

BUFFALO 

TWIN  OTTER 

78 

©  @ 

c2  -  -  <CL  +  CD  > 
a  a  o 

Cruise 

Cruise 

=  -  (5.2  +  .036)  =  -  5.24 

=  -  (5.2  +  .047)  =  -  5.25 

Slow  Fit 

Slow  -Fit 

=  -  (5.2  +  .057)  =  -  5.26 

=  -  (5.2  +  .065)  =  -  5.27 

Lndg  Appr 

Lndg  Appr 

=  -  (5.2  +  .127)  =  -  5.33 

=  -  (5.2  +  .155)  =  -  5.36 

79 

©  © 

Cz  =  -  2  a  V / 

q 

=  -  2  (3.5) (1.12)  =  -  7.83 

=  -  2  (3.5) (.915)  =  -  6.40 

80 

(p"N  fjp 

Cz.  ~  aT  d*  S  nlT 

e  (“f  ©(b) 

-  0.5)  (.54)  (233)  (1.0) 

=  .465 

=  (3.5)  (.54)  <i£2.)  (1.0) 

420  *>450 

81 

Cm  mV~  <h“hn> 
a  a  n 

=  (5.2) (-.15)=  -  .78 

=  (5.2) (-.15)  =  -  .78 

(for  normal  CG) 

(for  normal  CG) 

C16 


NO. 

BUFFALO 

TWIN  OTTER 

82 

S'© 

C  =  C  — 

m*  z .  e 

©  XD 

«  (-1.33)  (-^-)=  -  6.05 

10.1 

=  (-1.60)  (f^-)  =  -  6.15 

6.5 

83 

©  S'© 

% =  %  r 

XD 

=  (-7.83)  (-^-)  =  -  35.6 

10.1 

=  (-6.40)  (— )  =  -  24.6 

6.5 

84 

c.  .  7 

s«  s°  \j) 

46 

-  (.465)  (—  )  =  2.12 

10.1 

=  (.450) (  —  )  ■  1.73 

6.5 

85 

(not  used) 
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1 

93 

C  A  =  C  — 

: " 

10 

=  (-  .233)  (— )  *  -  .024 

96  , 

1 

1 

1 

=  (-.317)  (|g)  =  -  .024 

94 

d\  /©'y  P 

%  *  /  2  ««y 

■  2(5.2)  (.43)  j  (i4.82-.186y)ydy 
(945) (96)  ^3 

(  ; 

/ 

=. 0000493[7. 41(472-333) 

,  1 

-.062  (473-333)]  =  .20 

,  1 

-  2(5-2)('401/316.5yay 

(420) (65)  }4  yY 

/312  -  142  \ 

—  nnnoo  1  -  1  —  *>q 

1 

.  1  2  ) 

95 

1  A 

1  . 

cn  =  cn  +  cn  ' 

6  *F  ^  fuselage 

e 

"ft  Yfi  5“ 

r  F  XD 

■  -  •362>«H>-  •“%»© 

<volfu Sr  7 

CnB  1,3  Sb  w 

fuse  ©"  "©  & 

*  "  1,3  (945) (96) (“Io)"  ”-065 

C_  =  .166  -  .065  =  .101 

p 

=  -  (-  .492)  (^~)  =  .193 

i 

(1300)  7 

■-1-3  <420,  <65,^  ■-  -°72 

c  =  .193  -  .072  =  .121 
n3 
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I 
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BUFFALO 


TWIN  OTTER 


Cn  "  Cn  +  Cn 
P  Pw  /-n  Pp 


Lo  ri  _  La  *, 
~  11  iXR  ] 


iXR  J 

\D 


—  Ti  5,21 
4  L1  “  5T7TttJ 

*  -  .207  CL 


yD  b 


XD 


-  (-  -055) (||)  *  .025 


f  1  .  hi] 
L1  iottJ 


-  .208  CT 


-  -  (-  .085)  .033 


Cruise 


Cruise 


Cn  ^  -  .207 ( .30)  +  .025 
P 

«  -  .037 


-  -  . 208 ( . 42)  +  .033  -  -  .054 


Slow  Fit 


Slow  Fit 


-  .207 ( .77)  +  .025  -  -  .134  =  -  .208(.78)  +  .033*-  .129 


Lnd?  Appr 


Lndg  Appr 


-  -  .207(1.49)  +  .025  -  -  .283  =  -  .208(1.65)+  .033  *  -.310 
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cn  ”  cn  +  cn 

r  rF  rw 

.  ..V®.  >e 

S&  - 

=  -  («■>  (.368)  *  -  .169 

M 

s.  ■  •  <r 

=  -  (^-5)  (.429)  -  -  .168 

Cruise 

Cruise 

Cn  -  ^7^  -  .169  ■  -  .171 

Cn  =  -  -  .166  =  -  .171 

r  4 

r  4 

Slow  Fit 

Slow  Fit 

.025  , 

-  -  ,  .169  ■  -  .175 

4 

•025 

“  "  .  -  .168  =  —  .174 

Lndg  Appr 

Lndg  Appr 

-  -  .169  =  -  .188 

4 

•093 

*  ■  ^  .168  ■  ~  .191 

98 

nx  y,  b~ 

r 

-  -  (-  .233)  (i|)  =  .107 

=  -  (-  .317)  (2g-5)  =  .124 

99 

C  -  0 

“6a 

99 
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BUFFALO 


TWIN  OTTER 


NO. 

PARAMETER 

_ 

CRUISE 

SLOW 

FLIGHT 

CRUISE 

SLOW 

FLIGHT 

LANDING 

APPROACH 

100 

q  S 

132 ,500 

52,000 

26,850 

28,750 

15,500 

7,260 

101 

q  s/u0 

331 

242 

174 

103 

88 

61 

102 

q  SC 

1.34x10s 

.5252x10^ 

.2715X106 

185,500 

100,800 

47,100 

103 

q  sc/uQ 

3350 

2440 

1760 

667 

572 

392 

104 

q  Sb 

12.75xl06 

5.0X106 

2.58X106 

1.855X106 

1.008xl0S 

.471x10® 

105 

q  Sb/UQ 

31,850 

23,200 

16,700 

6,670 

5,720 

3,920 

106 

C/200 

.0126 

.0235 

.0328 

.0117 

.0185 

.027 

107 

b/2U0 

.120 

.223 

.312 

.1170 

.185 

.270 

108 

m 

1240 

1240 

1240 

372 

372 

372 

109 

*x 

273,000 

273,000 

273,000 

24,300 

24,300 

24,300 

110 

•  • 

I 

y 

215,000 

215,000 

215,000 

22,000 

22,000 

22,000 

111 

447,000 

447,000 

447,000 

41,000 

41,000 

41,000 

In  this  appendix,  transient  responses  (or  time  histories) 
of  aircraft  motions  in  response  to  step  control  inputs  are 
presented.  Equations  48-53  of  Section  V  and  the  numerical 
values  developed  elsewhere  in  this  report  are  used. 

Figure  D1  shows  the  response  in  pitch  rate  0,  pitch  angle 
A0,  angle  of  attack  a,  altitude  Ah,  and  forward  speed  (AU/UQ) 
resulting  from  a  step  elevator  input  A6e  for  the  “Buffalo" 
aircraft  in  the  cruise  flight  condition. 

Figure  D2  shows  the  response  (again  for  the  Buffalo  in 
cruise  flight)  in  the  same  parameters  resulting  from  a  step 
change  in  thrust  AT. 

Figures  D3  and  D4  present  Buffalo  A6e  and  AT  responses 
respectively  with  the  aircraft  at  the  slow  flight  condition. 

Figures  D5  and  D6  present  the  same  responses  for  the  landirg 
approach  condition.  In  Figures  D5-D6  (and  in  Figures  D17-D18) , 
altitude  (Ah)  in  plotted  as  the  change  from  the  nominal  altitude 
resulting  from  a  steady  descent  at  -7.5°  flight  path  angle. 

Figures  D7-D12  present  lateral  responses  for  the  Buffalo. 
Figure  D7  shows  the  response  in  sideslip  angle  0,  roll  rate  <J>, 
roll  angle  4>,  yaw  rate  $,  and  yaw  angle  4*  resulting  from  a  1° 
step  aileron  input  6  at  the  cruise  flight  condition. 


Figure  D8  shows  the  lateral  response  in  cruise  flight  due 
to  a  1°  step  rudder  input  Sr. 

Figures  D9  and  DIO  present  6a  and  6  responses  respectively 
for  the  Buffalo  at  the  slow  flight  condition. 

Figures  Dll  and  D12  present  the  same  information  for  the 
landing  approach  configuration. 

Finally,  Figures  D13-D24  present  the  same  information  for 
the  Twin  Otter  aircraft  as  is  given  in  Figures  D1-D12  for  the 
Buffalo . 

Time  constants,  natural  frequencies  and  damping  ratios 
associated  with  the  transient  responses  shown  in  figures 
D1-D24  are  presented  in  Table  Dl.  These  were  generated  as 
part  of  the  transient  response  computer  program. 
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TABLE  D1 

Summary  of  Characteristic  Response  Parameters 
For  "Buffalo"  and  "Twin  Otter"  Aircraft 


BUFFALO 


CRUISE  flight  APPROACH 


TWIN  OTTER 


CRUISE 


APPROACH 


SHORT  PERIOD  OSCILLATION 


in  _  RAD 

n  SEC 

2.93 

1.98 

1.42 

3.14 

2.46 

1.69 

c 

.794 

.855 

.856 

.710 

.780 

.783 

1/^U)n-SEC 

.43 

.59 

.82 

.45 

.52 

.76 

LONG  PERIOD  OSCILLATION 

mu 

.084 

.147 

.205 

.132 

.198 

.289 

.166 

.108 

.082 

.140 

.101 

.069 

l/£o>n-SEC 

72 

63 

59.8 

54 

50.3 

50 

LATERAL  OSCILLATION 

m  RAD 

wn  "  SEC 

1.78 

1.26 

1.09 

2.46 

1.95 

1.66 

.162 

.169 

.193 

.202 

.254 

.360 

1/Cwn-SEC 

3.48 

4.69 

4.77 

2.02 

2.02 

1.67 

SPIRAL  MODE  (MINUS  SIGN  INDICATES  DIVERGENCE) 


75.7 

-379 

-78.5 

788 

ROLL  SUBSIDENCE 


xR  -  SEC 

.328 

_ i 

.446 

.650 

.185 


-48.5  -21.8 


221  .376 


Oi-h0) ft/10 


FIG. D- 7  Response  To  1°  Step  Aileron  Input  (Buffalo, Cruise) 
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FIG.D-9  Response  To  1°  Step  Aileron  Input  (Buffalo, Slow  Flicht) 
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FIG.D-10  Response  To  1°  Step  Rudder  Input  (Buffalo, Slow  Flight) 
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FIG.D-16  Response  To  10%  Thrust  Input  (Twin  Otter,  Slow  Flight) 
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FIG.D-18  Response  To  10%  Thrust  Input  (Twin  Otter,  Approach) 
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TIME -sec 

FIG.D-20  Response  To  1°  Step  Rudder  Input  (Twin  Otter , Cruise) 
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